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ABSTRACT 



Expanding upon the IRAC survey from lDahm and Hillenbrandl (120071 ). we de- 
scribe Spitzer IRAC and MIPS observations of the populous, 5 Myr-old open clus- 
ter NGC 2362. We analyze the mid-IR colors of cluster members and compared 
their spectral energy distributions to star-|-circumstellar disk models to constrain 
the disk morphologies and evolutionary states. Early/intermediate-type con- 
firmed/candidate cluster members either have photospheric mid-IR emission or 
weak, optically-thin infrared excess emission at A > 24 /xm consistent with debris 
disks. Few late-type, solar /subsolar-mass stars have primordial disks. The disk 
population around late-type stars is dominated by disks with inner holes (canoni- 
cal 'transition disks') and 'homologously depleted' disks. Both types of disks rep- 
resent an intermediate stage between primordial disks and debris disks. Thus, in 
agreement with previous results, we find that multiple paths for the primordial- 
to-debris disk transition exist. Because these 'evolved primordial disks' greatly 
outnumber primordial disks, our results undermine standard arguments in favor 
of a < 10^ year timescale for the transition based on data from Taurus-Auriga. 
Because the typical transition timescale is far longer than 10^ yr, these data 
also appear to rule out standard UV photoevaporation scenarios as the primary 
mechanism to explain the transition. Combining our data with other Spitzer 
surveys, we investigate the evolution of debris disks around high/intermediate- 
mass stars and inves tigate timescales for giant planet formation. Consistent with 



Currie et al.l (j2008al ). the luminosity of 24 /im emission in debris disks due to 



planet formation peaks at ^ 10-20 Myr. If the gas and dust in disks evolve 
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on similar timescales, the formation timescale for gas giant planets surround- 
ing early-type, high/intermediate- mass (> 1.4 Mq) stars is likely 1-5 Myr. Most 
solar /subsolar- mass stars detected by Spitzer have SEDs that indicate their disks 
may be actively leaving the primordial disk phase. Thus, gas giant planet for- 
mation may also occur by ~ 5 Myr around solar /subsolar-mass stars as well. 

Subject headings: stars: pre-main-sequence-planetary systems: formation — 
planetary systems: protoplanetary disks - stars: circumstellar matter - Open 
Clusters and Associations: Individual: NGC 2362 



Facility: Spitzer Space Telescope 



Introduction 



Primordial circumstellar disks of gas and small dust grains s urround the vast major- 
ity of ~ 1 Myr-old stars an d are the building blocks of planets (iLada and Adamd Il992 



Kenyon and Hartmannlll995l ). Observable as emission in excess of the stellar photosphere 
at infrared wavelengths, primordial disks accrete gas onto the host star and have strong, 
optically-thick emission comparable in luminosity to the stellar photosphere (L^i/L^, ~ 0.1- 
1). Small dust grains in primordial disks comprise the seed material for the cores of planets; 
the gaseous envelopes of gas giant and ice giant planets form from the circumstellar gas in 
these disks. 

Primordial disks disappear on ~ 3-10 Myr timescales. In clusters with ages of ~ 3 
Myr, many disks have much less emission than typic al primordial disks at wavelengths < 
10 /im (e.g. iLada et al.ll2006l : iHernandez et al.ll2007ar). Pr imordial disks disappear around 
nearly all stars by ~ 10-15 Myr (e.g. Currie et ahlboOTa ). Disks surrounding > 10 Myr 



old s t ars are typically optically-thin c?e6ns disks w it h weaker infrared ern i ssion (IRieke et al. 



2005 


; 


2008 





Chen et al. 


2005; 


Currie et al. 


2007a. 


2008 Jbl: 1 



Plavchan et al. 



gas (e.g. 



Dahm 



20081 



and little evidence for substantial reservoirs of circumstellar 
Because dust from debris disks can be re- 
moved by stellar radiation on very short (5, 0.1 My r) timescales (e.g 



2005al : ICurrie et all l2007d ) 



Burns et al. 1979 



Takeuchi and Artymowica l200ll : iPlavchan et al.ll2005l ). the dust requires an active replen- 
ishment source from collisio ns between larger objects. Thus, debris disks are signposts fo r 
active planet formation (e.g. iBackman and Parescdll993l : iKenyon and Bromleyll2004l . |2008| ). 



While debris disks may trace active terrestrial and icy planet formation, they do not 
trace the final stages of active gas giant planet formation, where massive cores accrete 
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large gaseous envelope s . Debris disk have very low masses of residual gas, << 1 Mm, (e.g. 



Zuckerman et al.l Il995l : IChen and KampI |2004J : iRoberge and Weinbergerl l2008l : [France et al. 



20071 ). This mass is too small to contribute to a gaseous envelope with a mass comparable 
to those found in Jupiter and Saturn. 

The timescale for the disappearance of primordial disks and the subsequent dominance of 
debris disks has important implications fo r gas giant planet fo rmation. In the core accretion 
theory for gas giant planet formation (e.g. iPoUack et al.lll996l ). > Mars-mass protoplanetary 
cores can accrete a gaseous envelope, which is in hydrostatic equilibrium. However, once the 
cores reach f« 5-10 in mass, they can rapidly accrete a much more massive envelope. 
If the nebular gas disperses before the critical core mass is reached, g as giants fail to form. 



Typi c al timescales for massive c ores to forn a range from ^ 1 Mvr (IRafikov 



2004 



Currie 



20051 : iKenvon and BromlevI 120091 ) to 10 Myr flPoUack et al.lll996l : iRice and Armitagell2003f ). 
Because debris disks are gas poor, the time when most primordial disks evolve into debris 
disks sets an empirical upper limit for the formation timescale of gas giant planets. 

Identifying how primordial disks evolve into debris disks may also have important im- 
plications for planet formation. Early work from IRAS and ground-based data (specifically 
Taurus- Auriga) suggested that ~ 1-10% of ~ 1 Myr-old stars have so-called 'transition disks': 
disks with levels of mid-infrared emission in between those of stars surrounde d by primordial 



disks and stars lacking a disk as determined from color-color diagrams (e .g. iHartigan et al. 
1990l : ISkrutskie et~aDll99ol : Isimon and Pratolll995l : IWolk and Waltei]ll996h . Later studies of 
'transition' disks began to associate the term with a specific disk morphology: disks with 
inner regions that are optically thin and thus de pleted or devoid of gas and dust but colder, 
outer regions with optically thick emission (e.g. ICalvet et al.ll2005l ). If all disks follow make 
this transition, then disks are cleared from the inside out with a transition time of ~ 10^- 
10^ years. Plausible physical processes for removing primordial disks must then reproduce 
this rapid transition tim escale. UV photoevaporation models a r e successful in clearing disks 
within ~ 0.01-0.1 Myr (jClarke. Gendrin. and Sotomayorll200ll : [Alexander et al.ll2006l ). 



However, recent spectral energy distribution (SED) analysis of some 2-3 Myr-old stars 
indicates tha t they may follow a different evolutionary sequence from primordial disks to 
debris disks (iLada et al.l l2006l : ICurrie and Kenyoru |2009| ) . Their near-to-mid IR fluxes are 
consistent with a 'homologous depletion' in disk emission at all wavelengths through 24 /im 
instead of an 'inside out' evolution, indicativ e of a reduced ma ss of small (< 100 fim - 1 
mm) dust grains at all stellocentric distances (IWood et al. II2002I ) that could result from the 
coagulation of dust into 1 mm - Ikm-sized planetesimals and depletion of solid material by 
accretion onto the star. If many disks evolve from primordial disks to debris disks along 
this alternate path, then the primordial-to-debris disk transition timescale estimated from 
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the frequency of canonical 'transition disks' will be underestimated. Furthermore, if most 
primordial disks last longer than ~ 1 Myr, the frequency of 'transition disks' around 1 Myr- 
old stars may be low either because the transition timescale is fast or because most disks 
simply haven't had enough time to start leaving the primordial disk phase. If the latter 
is the correct interpretation, then the transition timescale for many disks could be much 
longer. Observations of disks at slightly greater, ~ 4-5 Myr ages can further test whether 
there are multiple pathways from primordial disks to debris disks and can better constrain 
the primordial-to-debris disk transition timescale. 

In this paper, we investigate the primordial-to-debris disk transition in the open cluster 
NGC 2362 in Cani s Maior (g^nnn = 7H8"'48 ^ (Jsooo = -24°57'00.0"). With an age of ~ 4-5 
Myr (jDahnj|2005bl : iMayne and Naylorll2008l ). NGC 2362 lies between the epoch when most 
stars have optically-thick primordial disks (~ 1 Myr) and the epoch where primordial disks 
are exceedingly rare (~ 10 Myr). While the cluster is very distant (~ 1.48 kpc), it is also 
massive (> 500 M©). Because of its mass, even nondetections of cluster stars with the MIPS 
24 /im channel can help con strain the source sp ectral energy distributions (SEDs) and thus 
the disk evolutionary state (ICurrie et al.ll2007bl ). 



This paper builds upon the study of NGC 2362 by iDahm and Hillenbrandl (120071 ) who 
focus on IRAC observations of 337 candidate/confirmed cluster members and find evidence 
for disks surrounding many late-type stars. Here we add MIPS 24 /xm data to extend the pre- 
vious Spitzer wavelength coverage of cluster members. We also rereduce and reanalyze IRAC 
data and expai id the membership sample so that we can include candidate cluster members 
not studied by iDahm and Hillenbrandl (120071 ). MIPS 24 /im data probes disk emission that 
can come from dust colder than the ice sublimation temperature (T ~ 170 K). Therefore, 
our knowldedge of the disk population now extends from the warm, inner terrestrial zone 
regions to cooler, outer regions where gas giant planets may form. The combination of MIPS 
and IRAC data allow us to better discriminate between disk evolutionary states. 

We organize our study as follows. In §2, we describe IRAC and MIPS photometry of 
the NGC 2362 field and ancillary data used to identify Spitzer counterparts to known and 
candidate cluster members. §3 focuses on the IRAC+MIPS colors of (candidate) cluster 
members, providing a first-order constraint on the mid-IR properties of disks in NGC 2362. 
To constrain the evolutionary state of MlPS-detected cluster members, we compare the 
observed SEDs to predictions for (evolved) primordial disks and debris disks in §4. §5 draws 
inferences about the timescale for gas giant planet formation and the primordial-to-debris 
disk transition from these data by comparing them with Spitzer results for other < 10-15 
MjTT-old clusters. We conclude with a brief summary of the paper, discuss our results in the 
context of other recent work, and discuss future work that can build on our results. 
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Data 



2.1. Spitzer Observations, Image Processing, and Photometry 

2.1.1. IR AC Data 



NGC 2362 was observed with the Infrared Array Camera (IRAC; iFazio et al.ll2004l ) on 
March 3, 2004 as a part of the Guaranteed Time Observations program (Program ID 37). 
IRAC data cover an area of ~ 0.3 square degrees centered near ^2000 ~ OT'^IQ'^OO*, ^2000 
~ -24°27'00". IRAC observations were taken in high-dynamic range mode with 0.4s and 
10. 4. s exposures. Sources were observed at two dither points yielding an average integration 
time/pixel of ~ 20.8s for the long exposures. 

The IRAC Basic Calibration Data (BCD) require additional post-BCD processing before 
mosaicing. Specifically, the BCD frames at [3.6] and [4.5] show strong evidence for column 
pulldown and 'striping' effects, which affect pixel flux densities by up to ~ 0.5 MJy/sr. The 
[5.8] and [8] channel BCD frames show evidence for cosmic-ray hits, which are especially 
numerous in the [5.8] filter. Finally, the automatically produced post-BCD mosaics at [3.6] 
and [4.5] contain background level mismatch of up to ~ 0.1 MJy/sr near bright stars, which 
causes a noticeable 'tiling' effect. 



We cleaned the IRAC data according to the methods des cribed bvlCurrie et al.l (l2008bl ). 
applying the array-dependent correction for point sources ( iQuijada et al.l 120041 ) and mit- 
igating image artifacts using modified versions of the muxstriping and column pulldown 
algorithms available as contributed software on the Spitzer Science Center website. Fi- 
nal image processing and m osaicing was performed using the MOPEX/APEX pipeline 
( iMakovoz and Marleaul |2005| ) . Prior to mosaicing, we applied the overlap corrector with 
bright source masking and 3a outlier clipping to match the background levels between BCD 
frames. For cosmic-ray rejection, we used the r-masks for each BCD frame supplemented by 
the mosaic outlier and dual outlier modules, which removed nearly all cosmic rays. Finally, 
both the IRAC and MIPS data were mosaiced together using a bicubic interpolation with 
outlier rejection. The final mosaics were inspected and cleared of image artifacts or pattern 
noise that could compromise the photometry. 

For the IRAC photometry, we per form source finding and a perture photometry using 
the IDLPHOT package and SExtractor (IBertin and Arnoutdl 19961 ) . For source identification, 
we determined the background rms level {(Jbkgd,rms) for each IRAC mosaic using SExtractor. 
Sources were identified as groups of pixels with fluxes > 3abkgd,rms in all filters using the 
find. pro routine. These sources lists were used as input for aperture photometry with the 
aper.pro routine. The source flux was computed in a 2 pixel (~ 2.44") and 3 pixel (~ 
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3.66") aperture radius with the local background calculated from a 4 pixel-thick annulus 
surrounding each source extending from 2-6 and 3-7 pixels. We multiplied the source fluxes 
by the appropriate aperture corrections given in the IRAC data handbook, version 3.cE|. We 
fine-tuned the aperture correction by first comparing photometry for bright, unsaturated 
(~ 10th-12th magnitude) stars derived above with photometry derived using a 10 pixel 
aperture and background annulus of 12 to 20 pixels, which requires no aperture correction. 
Additionally, we then compared our results with photometry for bright stars measured using 
a 5 pixel aperture with a 12-20 pixel background annulus. 

Photometry from the 2 and 3 pixel apertures show excellent agreement through ~ 13th 
(14th) magnitude for the [5.8] and [8] ([3.6] and [4.5]) channels with dispersions of << 
0.025 magnitudes. The dispersion increases to ~ 0.1-0.2 mags beyond these limits. The 
measured pixel area of > 13th magnitude stars is enclosed by the 2-pixel aperture radius. 
Using the 2-pixel aperture may introduce ~ 1-2% uncertainties in the spatial variation of the 
aperture correction but results in a higher signal-to-noise, which is important for obtaining 
accurate photometry of stars whose fluxes are < I0abkgd,rms- Furthermore, NGC 2362 is a 
massive cluster projected against the densely populated galactic plane. Source crowding in 
the sensitive [3.6] and [4.5] filters can be severe. Several stars in the cluster center appear 
partially blended with the three-pixel aperture used. Therefore, we chose photometry from 
the 2-pixel aperture in all filters for all stars. The final IRAC catalogs were trimmed of 
sources lying within 5 pixels of the image edges. We combine the IRAC catalogs with the 
2MASS JHKg catalog using a 2" matching radius, resulting in 6,916, 7,075, 4,875, and 3,855 
matches in the [3.6], [4.5], [5.8], and [8] channels 

The photometric errors are calculated from the Poisson error in the source counts, the 
read noise, the Poisson error in the background level, and the uncertainty in the backgrounc^l. 
Figure [T] shows the number counts of IRAC-detected sources with 2MASS counterparts as 
a function of magnitude. The source counts peak at [3.6] ~ 15.5, [4.5] ~ 15.5, [5.8] ~ 15.0, 
and [8] ~ 14.25-14.5. IRAC approaches the lOa limits at [3.6] - 15.75, [4.5] - 15.75, [5.8] - 
13.5, and [8] ~ 13.25. The 5a limits for [5.8] and [8] are [5.8] ~ 14.8 and [8] ~ 14.4. Source 
counts for [3.6] and [4.5] peak before the lOa limit is reached because of incompleteness in 
the 2MASS catalog. 



^http:/ /ssc.spitzer.caltech.cdu/irac/dh/iracdatahandbook3.0.pdf 

^See https://lists.ipac.caltech.edu/pipcrmail/irac-ig/2007-February/000037.html 
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2.1.2. MIPS 24 lim Data 
NGC 2362 was imaged with the Multiband Imaging Photometer for Spitzer (MIPS; 



Rieke et al.ll2004f ) on March 17, 2004 also as a part of the GTO program (ID 58). MIPS data 
cover ~ 0.5 square degrees on the sky with boundaries of 02000 ~ 07^17'^45'' to 07^20™'00'^ 
and ^2000 ~ -25''26/00" to -24°33/00". MIPS 24 /im observations were taken in scan mode 
with an integration time per pixel of 80s. Background cirrus levels at 24 /xm are low (^ 19 
MJy/sr) and vary by less than ~ 0.1 MJy/sr for most of the field. Two regions centered on 
^2000 ~ 07'^18'"8^ ^2000 ~ -25°16/00" and 02000 ~ 07'^19'"30^ ^2000 ~ -24°51/30" have slightly 
higher background levels of ~ 19.5 MJy/sr. The MIPS BCD data are largly free of artifacts. 

Because the MIPS point-spread function is well characterized, we perform photom- 
etry on the 24/xm MIPS data using pixel response function (PRF) fitting with APEX 



flMakovoz and Khan] l2005h . We detect 1,761 sources at [24] with S/N > 3. Using a 2" 
matching radius, we identify 447 MIPS 24 /im sources with 2MASS counterparts. Figure [2] 
shows that the source counts peak at [24] ~ 10.5 (~ 460 /xJy); the (5(t) lOa limits are [24] ~ 
(10.5) 9.7fl Table [Hhsts the full photometric catalog of 8,642 sources on the IRAC/MIPS 
fields with 2MASS counterparts and a detection in at least one IRAC and/or MIPS channel. 
Figure [3] shows a three-color image of NGC 2362 made with the 4.5 /im and 8 /im IRAC 
filters (blue and green) and MIPS 24 /xm filter (red). Most sources discussed in this paper lie 
within ~ 7/ of r CMa, near the center of the mosaic and interior to the region with slightly 
higher 24 /im background (righthand side of the figure). 



2.2. Ancillary Data: Identifying Cluster Members 

In the absence of proper motion data, we establish cluster membership by identifying 
stars with indicators of youth. Most intermediate and late -type young stars are chromo - 



spherically active and thus are luminous X-ray sources (e.g. iPreibisch and Feigelsorul2005l ). 
Chromospheric activity as well as circumstellar gas accretion onto a star can result in 
Hp, emission, a common feature in the op t ical spectra of star s with ages < 10 Myr (e.g. 



Kenvon and Hartmann 


1995; 


Dahm 


2005a: 


Currie et al. 


2007c) 



also identified as having strong Li 6708 A absorption. Finally, optical color-magnitude di- 



Lvra et al. 


2006; 


Currie et al. 


2009b) 



In this section, we describe ancillary data for NGC 2362 that identifies confirmed and 



•^The uncertainties quoted here do not include the zero-point uncertainty, which is « 4%. 
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cand idate cluster members. We dra w our primary list of confirmed and candidate members 
from IPahm and Hillenbrand! ( 20071 ) identified by Li 6708 A absorption, X-ray activity, Hq, 
emission and optical photometry. To this list, we add a seco nd list of candidate members 
identified by deep optical photometry from llrwin et al.l (120081 ). 



2.2.1. Optical Photometry and Membership List from lDahm and Hillenbrand 12002) 



Dahm and Hillenbrandl (120071 ) list 232 members from a sample of stars with Hq emission, 
Li 6707 A absorption, and X-ray activity. About 20 % of these stars show evidence of youth 
from all three indicators, ~ 36 % have at le ast two indicators of youth, a nd nearly half (~ 
44 %) have just one indicator. To this list, iDahm and Hillenbrandl (120071 ) add 105 sources 
whose optical colors/magnitudes are consistent with colors of ~ 5 Myr-old stars and are thus 
'candidate' members. While none of these stars have detected H^ emission, Li absorption, or 
X-ray activity, many of them are too faint (V > 21) to be studied with spectroscopy or X-ray 
surveys. Other stars are too early in spectral type for the nondetections of Hq, emission, Li 
absorption, or X-ray activity to rule out membership. 

Using a 2" matching radius, we identify all cluster members in at least one IRAC channel 
and 39 cluster members with MIPS. Increasing the radius to 3" does not yield any new MIPS 
matches and increases the likelihood of contamination by background galaxies. Thus, about 
12 % of the known cluster members are detected by MIPS: analysis of the mid-IR properties of 
cluster members (§3) is limited by many upper hmits . Table [2] lists the comb i ned op tical /IR 
(IRAC and MIPS) photometry of members from the iDahm and Hillenbrandl (120071 ) catalog, 
using IRAC photometry from this paper. 



While the membership list from lDahm and Hillenbrandl (l2007l ) is robust, it is spatially 
limited and does not include the lowest-mass stars. The grism spectra used to identify Hq, 
emitters are restricted to stars within a ll'xll' region surrounding the cluster center (near 
r CMa; a = 07^ 18™ 40.4^ 6 = -24° 33' 31.3"). Chandra/ACIS observations which detect 
identify X-ray a ctive n i embers are limited to a 16.9' x 16.9' area surrounding the cluster 
center. Finally, iDahml (l2005bl ) identifies candidate cluster members from photometry for 
stars with V < 21 within 7' of the cluster center. 



According to iDahml (l2005bl ). cluster members with H^ emission. X-ray activity, and/or 
Li absorption are very centrally concentrated. The surface density of members drops by a 
factor of two at a cluster radius of ~ 2' and m erges with backgrou nd by ~ 6 '-7'. How- 
ever, many massive clusters (e.g. h and x Persei; ICurrie et al.ll2007al ) contain a low-density 
halo/distributed population of young stars located many core radii away from the cluster 
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centers. The halo p opulation may contain as much mass as the bound cluster population 
dCurrie et al]l2009bh . 



2.2.2. Optical Photometry and Expanded Li st of Candidate Members from Urwin et al. 



To identify candidates for fainter members of the cluster and members of a surrounding 



halo-p opulation, we rely on the deep optic al photometry of NGC 2362 from llrwin et al. 



( 120081 ) . As a part of the MONITOR program JAigrain et al.lboOTi ) , llrwin et all (|2008|) imaged 
a 36' X 36' f ield centered on r CMa with the MOSAIC Imager on the 4m Blanco telescope at 
CTIO. The llrwin et al.l catalog includes ~ 56,000 stars with 5a detections at a magnitude 
limit of I ~ 23.6. 



Irwin et al.l (120081 ) identify ~ 1,826 stars brighter than V ~ 25 whose positions on 



the V/V-I color-magnitude diagram are con sistent with cluster membership. About 1,465 
of these stars lie outside the iDahnj (j2005bl ) survey limits, greater than 7' away from the 
NGC 2362 cluster center. All of these stars are fainter than V ~ 16.^. Using a simulated 



catalo g of field objects from the Besacon Galactic models (IRobin et al.ll2003l ). llrwin et al. 
(120081 ) estimate a contamination level of ~ 65 %, which implies that NGC 2362 contains a 
population of ~ 600-650 faint stars, ~ 500 of which lie in a halo population surrounding the 



cluster. Thus, the candidate member list from llrwin et al.l (120081 ) may identify new members. 



Because the llrwin et al.l (120081 ) catalog may include cluster members, combining their 
source list with the MIPS detections may reve al additional cluste r stars with mid-IR excess 
emission from disks. In total, 1,1 43 stars from Irwin et al.l (2008 ) have IRAC and/or MIPS 



counterparts and are not in the iDahm and Hillenbrandl (120071 ) membership sample. Nearly 
all (1,120) of these stars are located beyond the bound cluster region (< 7' from the cluster 
core). Table[3]lists the optical, 2MASS, and Spitzer photometry of these candidate members. 



^Brighter stars are saturated. 
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3. 2MASS/IRAC and MIPS Colors of NGC 2362 Cluster Members and 

Candidate Members 

To identify disk-bearing stars in NGC 2362, we analyze the mid-IR colors of confirmed 
and new candidate cluster members §. We use IRAC colors to identify stars with warm 
dust emission from the inner disk. Then we analyze MIPS colors to identify stars with 
colder emission. Comparing the IRAC and MIPS colors provides a first-order investigation 
of the NGC 2362 disk popul ation. Comparisons between our photometry and that from 



Dahm and Hillenbrand! (120071 ) are discussed in the Appendix. 



3.1. IRAC Colors 



Figures m and O show the distrib utions of K..-[5.8] (left p a nels) and Ks-[8] (right panels) 
colors for NGC 2362 stars f r om th e iDahm and Hillenbrand! (120071 ) sample (Figure Hj) and 
candidates from llrwin et al.l (120081 . Figure [5]) . To show how these colors vary with stellar 
properties, we plot them against the 2MASS J-K^ color. For cluster stars with similar ex- 
tinction, the observed J-K^ serves as a good tracer of spectral type. Late-type K and M stars 
have red J-K^ colors (> 0.6); early-type B and A stars have J- K.^ ~ 0. The positions of these 
stars in optical color- magnitude diagrams define a tight locus ( IDahm and Hillenbrandll2007l : 
Moitinho et al.ll200ll ). which suggests that the cluster lacks a large age spread. Therefore, 
the J-Ks color also correlates well with stellar mass. 

To compare the observed Ks-[5.8, 8] colors to predicted photospheric colors, we overplot 
the locus of photospheric colors (dotted lines) from the STAR-PET interactive tool available 
from the Spitzer Science Center websit^. All sources with photometric uncertainties of 
< 0.5 mags are shown. Sources with (without) MIPS detections are shown as black dots 
(triangles). 

In both panels of Figure HI the distribution of the IRAC colors follows the trends ex- 
pected for a young cluster. Most sources lie within ~ 0.2 mags of the predicted photospheric 
colors, which range from Ks-[8] ~ for early-type stars (J-K.. ~ 0) to ~ 0.5^0.6 for la t e-type 
stars (J-K^ ~ 0.8). Very few stars have K^-[5.8, 8] < 0. The IPahm and Hillenbrandl (l2007h 



^No source has a galaxy contamination flag in the 2MASS Point Source Catalog, indicating that its 
emission overlaps with that of an extended source. 

^ STAR-PET computes the 2MASS- S pitzer colors based on the Kurucz-Lejuene stellar atmosphere models 



Kurucz 


1993 


Leieune et al. 


1997 


)• 


Rieke et al. ( 



Rieke et al. I (|2008h shows that differences in IRAC and MIPS colors 



between synthetic and empirical A type and solar-type photospheres are < 0.05 mags. 
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sample lacks high-mass stars (J-K^ < 0.4) with red IRAC colors that are indicative of excess 
emission from warm dust. In contrast, many late-type, lower-mass stars (J-K^ > 0.6 - 
< 1.4 Mq at 5 Myr (IBaraffe et al.lll998l )) show evidence for IR excess emission. 



Most stars with clear IRAC excesses have MIPS detections. The MlPS-detected pop- 
ulation (black dots) has a narrower range in J-K^ (~ 0.8-1.4) than the IRAC-detected 
population. The IRAC colors of MlPS-detected stars also show evidence for a wavelength- 
dependent luminosity of disk emission. MlPS-detected stars characteristically have Ks-[5.8] 
~ 1. The Ks-[8] colors for most of these stars are slightly redder (Ks-[8] ~ 1.25-1.75). 



The distribution of IRAC colors for the llrwin et al.l (120081 ) candidates shows similar 
trends (Fig u re [5l) . The dispersions in K5-[5.8] and Ks-[8] colors are ~ 0.2 mags. The 



Irwin et al 



sample shows evidence f or a clear IR-excess populatio n (Ks-[5.8, 8 



> 



0.5-0.75), which is smaller than that for the iDahm and Hillenbrand! (120071 ) sample. MlPS- 
detected stars in this sample have Ks-[5.8] ~ 1 and Ks-[8] ~ 1-2. 



3.2. K5-[24] Colors and Upper Limits 



Figure M shows t he J vs. K..-[24] and J-K .^ vs. K .^- [24] diagrams for confirmed / c andidate 
cluste r members from lDahm and Hillenbrandl (120071 . top panels) and candi dates fromllrwin et al, 
(120081 , bottom panels) . Exti nction to the cluster is low (avg. E(B-V) ~ 0.1: lBalona and Laney 



1996 



Mathis 



Moitinho et al.l l200ll ) and affects optical fluxes far more than infrared fluxes (e.g. 
1990h. Therefore, potential 24 fxm excess sources are not affected by sm all dis- 



persions in the reddening of cluster stars. While the iDahm and Hillenbrand! (120071) sample 



contai ns some stars with Ks-[24] colors consistent with stellar photospheres, the llrwin et al. 
( 120081 ) sample lacks photospheric sources because it contains only faint stars. 



The Ks-[24] colors of confirmed/candidate members indicate that NGC 2362 harbors 
many stars with 24 fim excess emission plausibly due to circumstellar disks. MIPS detects 
13/22 of the known cluster members brighter than J=ll (Figure El top-left panel). These 
bright sources have J-K^ ^ 0, indicative of early spectral-type stars (Figure [6l right panels). 
Thus, the MIPS observations probably detect photospheric emission from many early-type 
(B-A) cluster members. The locus of sources with 24 /xm detections and upper limits in 
Figure E] (left panels) shows that MIPS detects both weak (Ks-[24] ~ 1-3) and strong (K^- 
[24] > 3-4) excesses around stars with J ~ 12-14.5 and only strong excesses around fainter 
stars. By comparing the locus of colors as a function of J magnitude and J-K^ color, we 
rule out the presence of strong 24 fim emission from cluster stars with J < 14, J-K^ < 0.6. 
For a 5 Myr-old cluster with negligible reddening, stars with J-Kg < 0.6 have spectral types 
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earlier than K3 and masses greater than ~ 1.4 ( iBaraffe et al.lll998l ) 



To summarize, the population of stars with 24 fim excesses reveals a wide range of mid- 
IR disk luminosities. Most bright, early-type stars have photospheric MIPS emission and 
thus may lack circumstellar dust. Early- type stars with MIPS excesses have weak excesses, 
^ 0.5-2 magnitudes. Later-type stars with MIPS detections typically have strong excesses. 
However, sensitivity limits preclude MIPS from detecting many late-type stars with weak 
excesses. 

The positions of several stars in Figure [6] call into question their me r abersh ip status. 
Specifically, IDs 28 and 166, listed as members in iDahm and Hillenbrand! (120071 ) based on 
X-ray activity, have J-K^ colors indicative of late-type stars (J-K^ ~ 0.6 and 0.8) but 
apparently have photospheric Ks-[24] colors. These stars have optical and 2MASS fluxes 
comparable to the earliest-type stars in the cluster and are brighter than other late-type 
cluster members by ~ 3 magnitudes. If these stars are true pre-main sequence members, the 
huge implied age spread is at odds with recent work show ing that NGC 2362 hkel y had a 
single, explosive burst of star formation ~ 4-5 Myr ago (e.g. lMayne and Nayloiil2008l ). These 
stars are probably either X-ray active foreground M stars, extremely high-mass, post-main 
sequence cluster members, or extragalactic sources. 



3.3. Combined IRAC and MIPS Colors of Cluster Members from 



Dahm and Hillenbrandl (120071 ) and Candidates from llrwin et al.l (120081 ) 



When combined with the IRAC data, the MIPS data reveal a diversity of mid-IR disk 
emission. Figure [7] shows the Ks-[5.8, 8] vs. Ks-[24] c olor- color diagrarn s for stars in the 
Dahm and Hillenbrandl (120071 ) sample (top panels) and llrwin et al.l (120081 ) sample (bottom 
panels). To compare the observed colors with colors typical of primordi al circumstellar disks, 
we overplot the colors from the median Taurus SED (iFurlan et al.ll2006l ) reddened by E(B-V) 
= 0.1. Most sources with MIPS upper limits (grey arrows) have K<,-[8] ~ -0.2-0.5. About 
18 stars have > 3.5 magnitude MIPS excesses. Comparing their K<.-[24] colors in Figure [7| 
to the distribution of Ks-[24] colors vs. J-K^ in Figure |4| shows that all stars with strong 
MIPS and IRAC excesses are late type, low-mass stars. 

Most sources with IRAC and MIPS excess emission have weaker emission than the 
Taurus sources with primordial disks. Because many of the 85 Taurus sources used to 
construct the median SED have evidence for inner holes/gaps (e.g. DM Tau, UX Tau), stars 
surrounded by full primordial disks should have Ks-[5.8, 8] colors comparable to or greater 
than those for the median Taurus SED. Only four NGC 2362 sources have IRAC excesses 



- 13 - 



exceeding that of the median Taurus SED. Most sources with excess emission have weaker 
emission than the median Taurus SED at both IRAC and MIPS wavelengths. Sources with 
IRAC and/or MIPS excesses weaker than the median Taurus SED, show a further diversity. 
Among the sources with 24 /xm excesses equal to or greater than the median Taurus SED, 
nearly half have IRAC fluxes in at least one filter that are reduced by a factor of ~ 2 from 
the median Taurus SED. 



4. Analysis of the NGC 2362 Disk Population 

The combined IRAC and MIPS photometry of confirmed/candidate cluster members 
reveals sources with a wide range of mid-IR excess emission consistent with a range of disk 
evolutionary states. In general, stars with blue J-K^ colors (< 0.4; early-type stars) either 
lack 24 nm excess emission or have weak 24 fim excess emission (K^ -[24] ~ 0.5-2). Stars 
with redder (J-K^ > 0.6; late-type stars) colors have strong 24 /im excess emission (K^- 
[24] ~ 3-7) presumably from optically-thick circumstellar disks. Furthermore, the combined 
IRAC and MIPS colors show that stars with strong 24 fim disk emission have a wide range 
of emission from inner disk regions probed by the IRAC 3.6-8 fim filters. The variety of 
mid-IR fluxes from disks surrounding NGC 2362 stars then warrants and enables detailed 
comparisons with predicted fluxes from disks spanning a range of evolutionary states. 

In this section, we explore the NGC 2362 disk population in more detail by using models 
to constrain the evolutionary state of disks surrounding early type, high/intermediate- mass 
stars and late type, low-m ass stars. By comparin g the observed SEDs of early-type stars to 



best-match SEDs from the lRobitaille et al.l (120061 ) radiative transfer models, we estimate the 



mass of emitting dust needed to reproduce the observed disk emission. We then compare 
this dust mass to dust masses for disks in several evolutionary states and derive the dust 
removal timescales. 

For late- type stars, our approach is twofold. First, we compare the source SEDs to the 
median Taurus SED and to predictions for warm debris disk emission. The median Taurus 
SED serves as an empirically -based fiducial primordial disk model. In the debris disk model 



( iKenyon and Bromley! l2004j ). debris emission from terrestrial planet formation is calculated 
at 0.7-2 AU from a 1 Mq star and tracked for ~ 100 Myr. The disk luminosity as a function 
of wavelength at 5 Myr is converted into a flux and scaled to the stellar photosphere. In 
these comparisons, we emphasize that no attempt is made to fit the observed SEDs. Rather, 
comparing the source SEDs to the median Taurus and debris disk models can constrain 
whether the SED is 'primordial disk like', 'debris disk like', or in between. For a s econd 



check on the disk evolutionary state, we fit the observed SEDs to the iRobitaille et al.l (120061 ) 
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models to determine the allowed ranges of dust masses, inner hole sizes, and disk flaring 
indices using the SED fitting tool of lRobitaille et al.l (120071 ) . From this analysis, we determine 
the relative frequency of primordial disks, debris disks, and disks in an intermediate stage 
surrounding late-type stars. 



4.1. Mid-IR Disk Emission from > 1.4 Mq-Mhss Stars 
4.1.1. SED Modeling of > I.4 Mq Stars 
To model the mid-IR emission from early-type stars with 24 /xm exc ess emission, we 



compa re th eir SEDs to the grid o f 200,000 radiative transfer models from iRobitaille et al. 



(l2006l ). The IRobitaille et al.l (120061 ) models predict the thermal emission of sources including 
contributions from stellar photospheres, infalling envelopes, bipolar cavities, and accreting 
circumstellar disks. The model parameters are varied for stars with a wide range of stellar 
temperatures. Dust masses are derived assuming a standard ISM size distribution (n(a) 
oc a~^'^). By selecting the best-matching SEDs, we estimate the masses of small (< 100 
fim-l mm) dust grains that are required to reproduce the observed level of emission. We 
then compare these masses to typical dust masses for debris disks and prim ordial disks. 

" M0 to 5x10-9 M0 (cf. IChen et al. 



2005 




Low et al. 


2005 


). TVT 


(cf. 


Andrews and Williamsl 



Low et al.l 120051). Tvpical d ust masses for primordial disks range from 10 ^-10 ^ Mq 



observed J-K^ colors. We model the SEDs of the four sources - IDs 31, 144, 303, and 314 - 
with K,-[24] = 0.5 - 2.5. 

Figure [8] shows SED model fits to the four early/intermediate- type star s with MIPS ex- 
cesses , the results of which are summarized in Table |H The best-matching IRobitaille et al. 
( 2006 ) models overpredict the observed fluxes from these sources at all mid-IR wavelength^, 
which indicates that the dust masses required to reproduce the observed disk emission 
through 24 fim are considerably less than the fitted masses. Three of the four sources - 
IDs 144, 303, and 314 - have significant 24 fim excesses above the model stellar photo- 
sphere (dotted line). The 24 fim flux from ID 31 lies just barely above the photosphere 
((K,s-[24])/ cr([24]) 4) so t his source has marginal evidence for a disk. The dust masses for 
the best-fit IRobitaille et al.l models are ~ 1-10 x 10"^° Mq, or ~ seven orders of magnitude 



less than typical masses of dust in primordial disks. Unless these sources have enormous 



'''Even though the "fits" are o bviously poor and clea rly overestimate the mass of dust required to repro- 
duce the observed emission, the Robitaille et al. ( 2006h grid lacks SEDs for disks with lower dust masses. 
Therefore, these are technically the "best-fit" models. 
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reservoirs of extremely cold dust, which would not contribute significantly to the observed 
24 fim emission, the total mass of dust in small grains is likely small compared to primordial 
dust masses as well. The extremely low required dust masses required imply that the disks 
are optically thin through 24 nm. These dust masses are inconsi stent with rnasses typical of 
primordial disks but are consistent with debris disk masses (e.g. IChen et al.ll2006l ). 



4.1.2. Dust Dynamics in Disks Surrounding > I.4 Mq Stars 



Next, we use dynamical arguments to show that, that dust surrounding high/intermediate- 
mass stars is most likely from debris disks (Table Hj). In optically-thin conditions, small dust 
grains are removed from the circumstellar environment by Poynting-Robertson drag and 
radiation pressure. The ti r nescale for dust remo val by Poynting-Robertson drag is (e.g. 
Backman and Parescdll993l : lAugereau et al.lll999l ): 



tp-ij(yr) ~ 705 



PsLa 



o 



< Qabs > L: 



(1) 



where tau is the distance from the central star, s is the radius of the grain, ps is the grain 
volume density and < Qabs > is the absorption coefficient. 

To derive the drag timescale for grains in the NGC 2362 sources, we consider typical 
parameters for tau, < Qabs >, and pg. These sources have negligible 8 pm excesses but 
clear 24 pm excesses. Thus, the typical grain temperature is likely ~ 100-200 K. For grains 
in radiative equilibrium, tau ~ (L*/Lo)°'^ (T/280 K)~^. Assuming an a ge of 5 Myr, earl y 
and intermediate- type stars in NGC 2362 have luminosities ~ 10-56 L© (jSiess et al.ll2000l ). 
Therefore, r ~ 6-40 AU for T ~ 100-200 K. We adopt r ~ 10-30 AU, roughly the mean 
distances for 2 00 K and 120 K dust in our sample. For grains with ~ 1 g cm~^ and 
< Qabs > ~ 1 JSurns et al.lll979l ). 



tp.Riji) ~ 7 X 10^ 



L' 



1pm' 'lOAU 



(2) 



For these massive stars, P-R drag removes small grains on timescales shorter than the 
~ 5 Myr age of sources in NGC 2362. For 10 pm grains orbiting stars with ~ 10-56 L©, 
the drag time is tp-R < 6.3 x 10^ yr. Grains with s < 10 am are probably por ous, with 



p. 



(e.g. 



Augereau et al.l 119991 ). Porous 



0.1 g cm~^ instead of our adopted ~ 1 g cm" 
grains have very short lifetimes, tp-R < 6.3 x 10^ yr. Thus, our derived upper limit for 
tp-R is conservative. 
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If the ratio of the force from radiation pressure to gravity (P) exce eds I/4J, radiation 



pressure can remove smal ler grains on faster, dynamical timescales (e.g. 
Strubbe and Chianej 120061 ): 



Burns et al. 1979 



3L^ < Qpr{a) > 



pr 



> 1/2. 



(3) 



WnGM^csp, 

This equation can be rearranged t o solve for the grain size below which radiatio n pressure 
can remove dust from the system (IBurns et al.l 1 19791 : iBackman and Parescd Il993l ) : 



rnax,iim 



< 1.14/im 



,Mq Igcrn' 



Ps 



^ Qpr ^ 



(4) 



Ps 



For stars listed in Table HI s^ax ranges from 5.8 to 23.8 /im, assuming < Qpr > ~ 1 and 
i 1 g cm~^. Using a more realistic dust grain porosity makes even more grains subject 
to blowout by radiation pressure. Therefore, grains responsible for producing the observed 
24 /im emission should be blown out unless they are replenished by collisions between larger 
bodies. 

As long as the gas density is low, radiation pressure can remove small grains on short 
timescales. When the grains are poorly coupled to the gas, radiation removes grains on 
timescales comparable to the local dynamical timescale, < 10^ yr. In gas disks with Vigas 
~ 1-10 M®, small (< 10-50 /xm) dust grains are well coupled to the gas. Even in these 
circumstances, they are pushed by radiati on pressure from k, 10 AU to hu ndreds of AU in 
~ 10^ years if the disk is optically thin (ITakeuchi and Artymowica I2OOII ). Disk emission 



originating at tens of AU then requires an active replenishment source and thus is likely 
second generation. 

The proper ties of disks surrounding ear ly-type NGC 2362 stars are highly suggestive 
of debris disks teackman and Parescel Il993l ). The small IR excesses in dusts imply that 
the disks have very low luminosities compared to their stars and have a low mass in dust. 
The rapid removal timescales for this dust, regardless of whether there exists residual gas, 
indicates that the dust is likely second generation. Disks with copious amounts of circum- 
stellar gas almost without exception have much stronger broadband near-to-mid infrared 
disk emission than that observed here. Thus, all known observed and inferred disk prop- 
erties fit standard debris disk characteristics and not primordial disk characteristics. More 
sensitive measures of circumstellar gas (e.g. optical echelle spectroscopy) would provide an 
even better test of our conclusion that these disks are debris disks. 



Strubbe and Chiane ( 20061 ) argue that grains more easily become unbound, on a dynamical timescale 



once (3 > 1/2, not 1. 
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4.2. SEDs of Solar and Subsolar-Mass Stars 



The Appendix shows th e atla s of SEDs for all late-ty pe (candidate) cluster members 
from IPahm and Hillenbrand! J2007h a nd llrwin et all J2008h with K^-[24] > 1. The two late- 
type stars in iDahm and Hillenbrand! (120071 ) with questionable membership status lack disk 
emission and are not shown. In each panel, we o yerplot the la error bars for IRAC and 
MIPS data. For the stars with spectral types from iDahml (j2005bl ). we overplot the SED of 
the stellar photosphere (solid line/squares) using a model specific to the spectral type. We 
calculate the range i n reddening and extinct i on fr om comparing the observed V-I colors to 
intrinsic colors from iKenyon and Hartmanru (119951 ). assuming a one-subclass uncertainty in 
the spectral type. All but five sources have a range of extinctions consistent with E(B-V) ~ 
0.1 (IMoitinho et al.ll200ll ). For these five sources, we use a value for reddening halfway in 
between the minimum and maximum Ay values allowed. For the other sta rs with spectral 
types , we assume A„ = 0.31, which is the mean extinction to the cluster (IMoitinho et al. 



20011 ). For the stars without known spectral types, we fit the source SEDs from I through 
K-band with stellar photosph eres to determine the best-fit spectral type, using colors from 
Kenyon and HartmannI (119951 ) and assuming a reddening of E(B-V) = 0.1 (A„ ~ 0.31). 



Because most Taurus sources lack clear J-band excess emission (IFurlan et al.ll2006l ). we 
scale the median Taurus SED (dotted line) by its J-band fiux, setting it equal to the stellar 
photosphere. The median Taurus SED clearly has the strongest disk emission. The warm 
debris disk models (dashed line/ crosses) have weaker disk emission at all mid-IR wavelengths, 
slightly in excess of the stellar photosphere. The nominal MIPS 24 fim detection limit in 
XFx space is overplotted as a grey dashed line (AF^ ~ 5.75 x 10~^^ W cm~^). 

To compare the inferred physical states of disks with more simple empirical classifi- 
cations, we measure the IRAC [3.6]-[8] fiux sl ope for eac h star. We divide the sample 
according to the three empirical divisions from iLada et al.l (120061 ). We label sources with 
a > -1.8 as those with "strong IRAC" emission, those with a = -2.56 - -1.8 as having 
"weak IRAC" emission , and those with a < -2.56 as those with "photospheric" emission. 
Dahm and Hillenbrand! (120071 ) use this fiux slope to infer the physical structure of the disk, 
specifically identifying sources with a > -1.8 as having primordial disks. Comparing the 
empirical disk classifications w ith evolutionary states d etermined from SED modeling allows 
us to compare our results with iDahm and Hillenbrandl 's and test the utility of using a single 
IRAC slope to probe disk states. 

Tab le |5] summarizes our resu l ts and is divided into three samples: Ml PS-detected sources 
from the lDahm and Hillenbrandl (120071 ). MlPS-detected sources from the llrwin et al.l (120081 ) 
sample, and sources la c king clear MIPS detections that have disk emission according to 
Dahm and Hillenbrandl (120071 ). We list basic sources parameters (spectral type, extinction) 
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as well as derived di sk parameters - disk mass, minimum inner disk radius, and flaring 
index {j3) - using the iRobitaille et al.l (120061 ) models which, along with cor nparisons to the 



Taur us SED and debris disk SED, determine the disk evolutionary state (e.g. IRobitaille et al. 



20071 ). For each source, we consider all models which satisfy the goodness-of- fit criterion 



~Xhest < ""-data rather than simply the best-fit model. The derived disk evolutionary states 
for each sample are discussed in the next two sections. For the simple comparisons to the 
Taurus and debris disk SEP, we assurn e a sin gle extinction. The extinction is treated as 
a free parameter in the IRobitaille et al.l (120061 ) models. However, the derived evolutionary 
states are insensitive to changes in the assumed extinction over the range of plausible values. 



4-2.1. SEDs of Confirmed/ Candidate Cluster Members with MIPS detections. 



The SEDs of MlPS-detected sources shown in the Appendix reveal disks with a variety 
of morphologies. Figure [H] shows comparisons between the Taurus and debris disk SEDs with 
four representative SEDs from the atlas illustrating the range of evolutionar y states listed in 



Table [5] Figure [TO] shows modeling results for one SED (ID 267) using the IRobitaille et al. 



( 120061 ) models. Below, we describe the characteristics for each evolutionary state and discuss 
their relative frequencies: 



1. Primordial Disks - Four stars in the Dahm and Hil. 



Ill, 137, 139, 202) and two stars in the llrwin et al 



enbrandl (120071 ) sample (IDs 



mm sample (ID-IR526 and 



ID-IR1046) likely harbor primordial disks. These sources typically have near-to-mid 
IR fluxes comparable to the median Taurus SED, though ID 202 appears to have a 
slight near-IR flux deficit. Only one of the six sources has an SED with mid-IR clearly 
stronger than the median Taurus SED (ID 137). Its stronger emiss ion may be due 



to a g reater level of disk flaring or inclination. Comparisons with the IRobitaille et al. 



(120061 ) models show that these disks probably do not have inner holes and most have 
masses of M^j^fc/Mo ~ 10~^-10~^'^ for a solar gas-to-dust ratio, comparable to both 
the Minimum Mass Solar Nebula an d the typical masses of primo rdial disks in Taurus 
as derived from submillimeter data (lAndrews and Williamsll2005l ). 



2. Homologously Depleted Disks - Ma n y sta. rs in both the iDahm and Hillenbrand 



( 120071 ) sample (12 stars) and llrwin et al.l (120081 ) sample (5 stars) show evidence for a 
disk that is significantly depleted in small dust grains through 24 /xm relative to typical 
primordial disks: a 'homologously depleted' disk. Like the median Taurus SED, the 
SED slopes for these sources typically follow a nearly power-law decline with wavelength 
through 24 /im. However, these sources typically show a flux deficit in the IRAC and 
MIPS bands relative to the median Taurus SED. Indeed, they are underluminous by 
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a factor of > 2-3 from 5.8 fim to 24 /im. Thus, they have a steeper power-law dechne 
of disk emission with wavelength compared to primordial disks: AF^ oc A"" where n = 
-7/4 to -5/4. These disks also lack inner holes but have dust masses in small grains that 
are ~ 100-1000 times less massive than their primordial disk counterparts The dust 
masses in small grains required to reproduce the observed SEDs are 100-1000 times 
less than those of their primordial disk counterparts. They also often have little flaring 
(with a flaring parameter, /3, of ~ 1-1.1), which is consistent with a significant amount 
of dust settling. These features are consistent with a disk that loses a significant 



amou nt of mass of small dust grains at all disk radii simultaneously (e.g. IWood et al 



2OO2I ). which could plausibly be due to the growth of dust to larger, 1 mm-1 km-sized 



bodies and/or accretion onto the star. 



3. Transition Disks - Five star s in the lDahm and Hillenbrandl (120071 ) sample and seven 
stars in the llrwin et al.l (120081 ) sample have SEDs consistent with the presence of inner 
holes/gaps indicative of so-called transition disks. Compared to the median Taurus 
SED, these sources have weaker IRAC fluxes but comparable MIPS 24 fim fluxes, a 
featu re typical of disks with inner holes/g aps suc h as DM Tau, H e n 3-6 00A, and TW 
Hya JCalvet et all I2OO5I : ILow et al.ll2005f ). The iRobitaille et"aD J2006h models that 
provide a good fit have minimum inner holes ranging in size between 2.6 AU and 7.1 
AU. Their dust masses range between ~ 10~^ Mq and 10""^ Mq. For a solar gas-to- 
dust ratio, the inferred disk masses range between ~ lO'^ M© and ~ 10"^ Mq. These 
values indicate that the disks are slightly depleted relative to primordial disks. One 
source (ID 194) has strong K-band and 24 fim excesses but weak IRAC excesses. Its 
SED may be suggestive of a disk gap instead of an inner hole. 



While most sources have SEDs that are straightforward to classify, several have charac- 
teristics that make them slightly more ambiguous. Several sources (IDs 41, 63, and 196) have 
photospheric emission shortwards of 5.8 /im or 8 fim; their 24 fim excess emission is much 
weaker than other late-type stars detected by MIPS (Ks-[24] ~ 0.75-2.5). In this sense, 
their disks are consistent with being debris disks. However, dust grain removal processes 
(e.g. Poynting-Robertson drag) likely operate on far longer (~ 1 Myr) timescales in disks 
surrounding these stars than disks surrounding high and intermediate-mass stars examined 
in §4.1.1. Therefore, that dust in disks surrounding these stars need not be be second genera- 
tion debris dust. Based on their dust masses and presence/absence of inner holes, we classify 
these sources as either homologously depleted disks or transition disks instead. Identifying 
these sources as having debris disks changes the disks statistics we discuss later but will not 
qualitatively change our results. 
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Most (15/21) of the iDahm and Hillenbrand (120071 ) sample with large MIPS excesses 
have measured Hq, equivalent widths. Because Hq, emission can be due to circumstellar gas 
accretion, its presence also provides some constraints on disk evolutionary stat e s. Us ing 
the results of iDahm and Hillenbrand! (120071 ) and the criteria of IWhite and Basril (120031 ) to 
identify accretors, two of the three sources with primordial disks and Hq, measurements show 
evidence for strong Hq, emission and thus are accreting. None of the sources with extremely 
low disk masses (IDs 41, 63, and 196) show this accretion signature. Stars surrounded by 
transition disks and homologously depleted disks include both accretors and those lacking 
evidence for accretion. Not all sources have optical spectra, so it is difficult to make definitive 
statements about the relative frequency of accretion for each disk state. However, the lack 
of evidence for accretion for the sources with the weakest disk emission is consistent with a 
depleted reservoir of circumstellar gas in these systems. The low frequency of accretion for 
homologously depleted disks (20%) and their low mean Hq, equivalent widths ( mean value 
is 9.2A) is expected if these systems have lower rates of accretion and are generally more 
evolved than their primordial disk counterparts. 



4.2.2. Disk-bearing Stars in NGC 2362 Lacking MIPS Detections 



Dahm and Hillenbrandl (120071 ) identify 14 stars with strong IRAC- excess emission and 33 
with weak IRAC-excess emission. MIPS detects 16 stars included in the Dahm and Hillenbrand 



(I2OO7I ) list, leaving 31 without detections. Of the sources lac king MIPS detections, three h ave 
strong IRAC excesses a nd are classified as 'primordial' by Dahm and Hillenbrand ( 2007, )R: 
Dahm and Hillenbrandl classify the remaining 28 sources as 'weak'. Disk candidates without 
MIPS detections 1) could have 24 /im fiuxes fainter than the MIPS detection limit, 2) could 
be close to a bright 24 /im source, or 3) might not be true disk-bearing stars. 

To infer the physical state of disks surrounding stars without MIPS detections, we first 
examine the processed MIPS mosaic image to identify any stars that are ~ 2-3 a above 
the background or are blended with another star. Next, we identify stars that are accreting 
and thus either have primordial or more evolved primordial disks (transition/homologously 
depleted). Finally, we use the MIPS upper limits to constrain the possible range of 24 /im 
fiuxes and thus constrain the disk evolutionary state. 



Dahm and Hillenbrandl (|2007t ) classify disks as 'primordial' or 'weak' according to their IRAC flux 



sl opes. Thus , these classifications are empirical classifications. The division between classes follows that 
of Lada et al. ( 20061 ). Usage of the term 'primordial' in this paper always refers to the inferred physical state 
of the disk from comparing its emission to disk models and the median Taurus SED. 
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The bottom section of Table [5] lists the results of our analysis and Figure [TT] shows SEDs 
of five additional primordial and evolved primordial disk candidates. One star, ID 100, is 
likely surrounded by a primordial disk given its strong IRAC excess emission. In the MIPS 
24 fim filter, it appears to be blended with a brighter cluster member. Four stars - IDs 
12, 27, 39, and 227 - have clear IRAC excess emission and are barely discernable on the 
MIPS mosaic at very high contrast. The signal-to-noise for these sources is likely ~ 2-3, 
which would exclude them from the MIPS source list. ID 12 shows evidence for accretion. 
The SEDs of these four stars lie below the lower quartile for the median Taurus SED, which 
indicates that they are probably most similar to homologously depleted disks. The remaining 
sources lack evidence for accretion, are not discernable on the MIPS mosaic, and have weak, 
< 1 magnitude IRAC excesses. These sources either not are true disk candidates, have very 
large inner holes, or harbor optically-thin disks, possibly debris disks. 



4-2.3. Summary of SED Analysis 



Analyzing the SEDs of late type, MlPS-detected NGC 2362 stars reveals disks with a 
wide range of evolutionary states. Primordial disks in N GC 2362 are very rare. Only 1 9% of 
MlPS-detected confirmed/candidate cluster stars fror n iDahm and Hillenbraridl (120071) have 
SEDs co nsistent with prina ordial disks. Combining the lDahm and Hillenbrand! (120071 ) sample 
with the llrwin et al.l (120081 ) sample of 14 MlPS-detected stars indicates that the true fraction 
of primordial disks in NGC 2362 is even lower, ~ 17% (6/35). 

M ost disks surrounding late-typ e cluster stars are in an advanced evolutionary state. 
In the bahm and Hillenbrandl J2007h sample, ~ 57% (12/21) of late type, MlPS-detected 
sources have SEDs consistent with disks wi th a reduce d mas s at all disk radii (homolo- 
gously depleted disks). If candidates from llrwin et al.l (120081 ) are included, the fraction 
of these disks is ~ 49% (17/35). Disks with inner holes/gaps ('transition disks') also 
comprise a large fraction of Ml PS-detected sources, representing ~ 24% (5/21) of MlPS- 
detected [DaLiiosijiillenbrai]^ sou^^^ and ~ 34% (12/34) of all confirmed/candidate clus- 
ter members detected by M IPS. In total, 81% (83%) of the MlPS-detecte d late-type stars in 
the iDahm and Hillenbrandl sample (IDahm and Hillenbrandl - Hlrwin et al.l samples combined) 
have disks in evolutionary states between primordial disks and debris disks. Any other disk- 
bearing stars lacking MIPS detections have 24 /zm emission lying below the detection limit. 
Compared to the stars with MIPS detections, disks around these stars have even weaker 
emission suggestive of substantial evolution. Therefore, disks around most late type 
stars in NGC 2362 are leaving the primordial disk phase. 



Our results stand in contrast to those from 



Dahm and Hillenbrand 



(l2007h who find that 
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primordial disks comprise ~ 34% of the total disk population. Our photometry identifies a 
comparable number of confirmed members with a > -1.8 (16) and between -1.8 and -2.56 
(30), so poorer photometry is unlikely to explain this discrepancy. However, many sources 
(e.g. ID 3) listed in Table E] and identified as 'primordial' disks by iDahm and Hillenbrand 
clearly have SEDs that do not look anything like canonical primordial disks found in Taurus 
and other young star-forming regions. In part, some of these descrepancies result from the 
adopted definition of a primordial disk. In this paper, we define a primordial disk to be a 
disk with an SED close to the median Taurus SED and a disk mass > 10^'^ Mq. Modeling 
shows that only five of the 16 stars with a > -1.8 have SEDs and disk masses typical of 
primordial disks. Using only a single parameter (the IRAC slope) then results in a larger 
number of primordial disks than identified from full SED modeling. 



5. Implied Constraints on The Evolution of Circumstellar Disks and Planet 

Formation 

The disappearance of primordial disks and the emergence of debris disks sets an empir- 
ical limit for the formation timescale of gas giant planets. The disappearance of IR excesses 
in disks implies the disappearance of the building blocks - solid material - for the cores of gas 
giant planets. Similarly, the lack of gas accretion implies a lack of circumstellar gas in inner 
disk regions, which corresponds to a lack of material to build the gaseous envelopes of gas 
giant planets. Therefore, in this section we use results from analysis presented in §4 to con- 
strain the evolutionary states of disks and to investigate the formation timescale of gas giant 
planets as a function of stellar mass (§5.1). We divide the sample into high/intermediate- 
mass stars {Mi, > 1.4 Mq) and solar /subsolar-mass stars (< 1.4 Mq). 

From the population of NGC 2362 stars surrounded exclusively by debris disks, we 
compare their 24 iim debris luminosity with that for > 5 Myr-old stars to investigate the 
evolution of debris emission and the connections this evolution has with the growth of icy 
planets (§5.2). In §5.3, we compare the more evolved primordial disk population to disks in 
other 1-10 Myr-old clusters to examine the morphology of and timescale for the primordial- 
to-debris disk transition. 

The results presented in this section are based on analyzing the MIPS- detected popula- 
tion of NGC 2362 confirmed/candidate members and other cluster stars with clear evidence 
for disks listed in Table [51 Most faint confirmed/candidate cluster stars are not detected by 
MIPS because of its bright 5a limit ([24] ~ 10.5). Therefore, completeness at 24 nm may 
affect our results. However, we show that any disk-bearing sources not detected by MIPS 
either harbor debris disks or homologously depleted disks with low masses. MIPS is sensitive 
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enough to detect levels of 24 fim emission typical of primordial disks surrounding faint clus- 
ter stars and emission typical of all but low- luminosity debris disks around brighter cluster 
stars. Thus, our ability to determine when primordial disks disappear is not undermined by 
sample incompleteness. Because different dust masses are implied for detections as opposed 
to non detections, the detection limit may serve as a discriminator between disk types. 

Because not all of the MlPS-detected members and candidate members have published 
spectral types, this section uses the J-K^ color as a tracer of spectral type and thus, assuming 
a 4-5 Myr age, stellar mass. While the lack of published spectral types for some sources adds 
uncertainty to our analysis, known properties of NGC 2362 stars show that this uncertainty 
is not significant given the coarse division into stars wi th masses greater than o r lesser 



than 1.4 Mq. NGC 2362 has a small dispersion in age (iMayne and Naylorl |2008| ) . which 
yields a reliable conversion from spectral type to stellar mass. The cluster also has a very 
low reddening (E(B-V) ~ 0.1), which translates into E(J-Ks) ~ 0.05-0.06. For early type, 
high/intermediate- mass stars (BAFG) to be misidentified as low-mass stars based on their 
observed J-K^ colors, they must have a reddening ~ 3-5 times the typical reddening of the 
cluster. There is little evidence for such a large spread in reddening. Cluster stars with 
known spectral types have J-K^ colors consistent with using J-K^ = 0.6 as the division 
between high/intermediate-mass stars and low-mass stars. 



5.1. Constraints on the Formation Timescale for Gas Giant Planets 



Assuming an age of ~ 4-5 Myr, NGC 2362 sta rs with masses > 1.4 Mq have G or earlier 
spectral types and J-K^ < 0.6 (IBaraffe et al.lll998l ) . Cluster stars with masses > 2 Mg^ have 
spect ral types earlier than mid F and J-K^ ^0.3 (IKenyon and Hartmannlll995l : ISiess et al. 
20001 ). From Figure [6], the Ks-[24] upper limits for stars with My, ~ 1.4-2 Mq are ~ 2.5; the 
upper limits for more massive, 2 Mq stars are ~ 1-1.5. In comparison the predicted Ks-[24] 
excesses for primordial disks and transition disks discussed in §4 are ~ 4-5.5 magnitudes 
over this spectral type range. Thus, any primordial disks or transition disks around cluster 
stars with masses > 1.4 Mq are easily detectable. Stars with 'homologously depleted' disks 
listed in Table [5] have Ks-[24] > 3 and are also easily detectable. 

Analysis from §4 shows that there are no primordial disks, transition disks, or homol- 
ogously depleted disks surrounding NGC 2362 stars plausibly in this mass range. MIPS- 
detected stars with M^ > 2 Mq have photospheric emission or < 2 magnitude excesses (cf. 
Figures [6] and [7]) . MlPS-detected stars with masses likely between 1.4 Mq and 2 Mq include 
only one star, ID 144, with nonphotospheric emission. High/intermediate- mass cluster stars 
with optical spectra do not show any evidence for Hq, emission indicative of gas accretion. 
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Based on SED modeling and dust removal timescales (§4), all disks surrounding high/intermediate- 
mass cluster stars are debris disks. Furthermore, even though Spitzer observations only trace 
thermal emission from small circumstellar dust grains not gas, the observable signatures of 
disks in evolutionary states with gas (primordial, transition, homologously depleted) are 
clearly lacking from any high/intermediate- mass star studied here. NGC 2362 then lacks 
any disks around high/intermediate- mass stars which have yet to reach the debris disk phase. 
Recent combined Spitzer and ground-based surveys of young clusters indicate that circ um- 
stellar gas should disperse on similar or faster timescales than dust (ILada et al. I booeh . If 
this is the case, our results suggest that gas giant planet formation likely occurs prior to 4^5 
Myr for > I.4 Mq stars. 

MlPS-detec ted late-type stars i n NGC 2362 hav e inferred masses between ^ 0.3 Mq 
and 1.4 Mq (cf. ISaraffe et aPllQQsl : Isiess et al.l[200ol : IPahm and Hillenbrandl l2007h . Most 
late-type stars detected by MIPS have SEDs consistent with more evolved primordial disks 
or debris disks. About 1/3 of MlPS-detected disks in the iDahm and Hillenbrandl (120071 ) 
sample are transition disks. Nearly half of all MlPS-detected stars are surrounded by disks 
that are depleted in mass at a range of disk radii (homologously depleted) . There is no clear 
connection between disks with this morphology and active planet formation. Only four stars 
have SEDs consistent with primordial disks. Given the MIPS detection limit for this sample, 
we can rule out the presence of any > 4.5 magnitude excesses typical of primordial disks and 
transition disks. 

Disks surrounding late-type stars in NGC 2362 that are not detected by MIPS must have 
Ks-[24] < 4-4.5. These disks are likely either debris disks or homologously depleted disks 
with levels of emission slightly below the MIPS detection limit. The disk population of solar 
and subsolar-mass stars is dominated by the large number of MIPS upper limits, making 
statistically robust conclusions about the bulk disk p opulation difficult. However, combined 
with the low frequency of accretion in NGC 2362 (IDahnj l2005bl : iDahm and Hillenbrand 
20071 ). these results demonstrate that most disks surrounding late- type stars have signifi- 
cantly evolved and that most lack evidence for circumstellar gas needed to form gas giant 
planets. 

While some primordial disks remain around < 1.4 Mq stars, the disk population for these 
stars is dominated by disks which have clearly undergone a significant amount of evolution 
(transition disks, homologously depleted disks). The morphologies of transition disks and 
homologously depleted disks are consistent with disks actively leaving the primordial disk 
phase for the debris disk phase. If these disks are actively making the primordial-to-debris 
disk transition, then gas giant planet formation also likely ceases by ~ 5 Myr for many, 
perhaps most, solar /subsolar-mass stars. 
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5.2. Debris disk evolution around High/Intermediate-Mass Stars 



Combining NGC 2362 data for high/intermediate- mass stars with data from otlier clus- 



ters probes tlie time evolution of debris emission from planet 
exara i ne the time evolution of 24 nm excess ern ission (see a lso 



20061 . [Hernandez et al.ll2006l . ICurrie et al.ll2008al . and IWvattll2008f ). focusing on 5-40 Myr- 



brmation. Specifically, we 



Rieke et al. 2005, Su et al. 



old early-type stars. We include all confirmed/candidate high/intermediate- mass members 
in NGC 2362 with MIPS detections. In addition to our data, we include sta rs from clusters 



with debris disks: O rion OBla and b (10 and 5 Myr; iHernandez et a 



Gautier et all 120071). the 3 Pic Moving Group (12 Myr; iRebuU et al 



■20061) . 7] Cha (8 Myr; 
2008h. h and y Persei 



(~ 13 Myr; ICurrie et al.ll2008af ). Sco-Cen, and NGC 2547 (38 Myr; iRieke et al.ll2005[ ). For 
the Sco-Cen, we include data for Upper Sc orpius (5 Myr), Up per Centaurus Lupus (16 Myr), 
and Lower Centaurus Crux (17 Myr) from lChen et al.l (120051 ). We limit our literature sample 
to BAF stars. 



Figur e [T^ shows the indiv idual amplitudes of 24 /im excesses as a function of time. In 
support of ICurrie et al.l (j2008al ). there is a rise in the level of debris e mission frorn ~ 5 M yr to 
10 Myr and a peak at ~ 10-20 Myr (see also ICurrie et al1l2008bl ). iRieke et al.l J2OO5I ) show 
that the emission at later ages (> 30-50 Myr) is consistent with a 1/t decline (dot-dashed 
line). The mean level of excess (dashed line) for < 20 Myr-old clusters clearly deviates from 
a t^^ decline. 

MIPS data for NGC 2362 s tars supports the tr end of debris disk evolution from 5 
Myr to 10 Myr first identified by ICurrie et al.l (l2008a] ). Three clusters/associations in Fig- 
ure [12] have ages of ~ 5 Myr: Orion OBlb, Upper Scorpius, and NGC 2362. All three 
clusters/associations have levels of excess that are weaker than those for ~ 10-15 Myr-old 
clusters, Orion Obla, the P Pic Moving Group, h and x Persei, and Upper Centaurus Lupus. 
Debris disks in 5 Myr-old clusters also have l ower 24 p,m luminos i ties than those in 25-40 
Myr-old clusters NGC 2232 and NGC 2547 JCurrie et al. l l2008bl : iRieke et al.ll2005l ). The 
conclu sion that debris disk s surrounding early-type, high-mass stars exhib it a 'rise' and 'fall ' 
trend (ICurrie et al.ll2008al ) instead of a monotonic decline in emission (e.g. iRieke et al.ll2005l ) 
is now very robust. 

The observed behavior of debris disk evolution agree s well with predictions from sophis- 
ticated debris disk models. iKenyon and Bromleyl (120081 ) investigate the evolution of debris 
emission from icy planet formation in disks surrounding > 1 Mq stars. For ~ 2 Mq stars 
with a disk mass of ~ 2-3 x a scaled Minimum Mass Solar Nebula, 24 /im debris emission 
rises rapidly from ~ 5 Myr to ~ 8 Myr after the start of the simulation which is sustained 
for ^20 My r. Spitzer observations indicate a ~ 3 Myr lifetim e for pr imordial disks around 



these stars (IHernandez et al.ll2007al : ICurrie and KenyonI 120091 ). The iKenyon and Bromley 
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models would then predict an observed peak in 24 /xm debris emission by ~ 8-11 Myr that 
lasts through ^ 30 Myr. These expectations agree with the trend shown in Figure | 



Within the context of the lKenyon and Bromleyl (120081 ) models, the trend in 2 4 fim debris 
emission is directly linked to the growth of icy planets. The lKenyon and Bromleyl simulations 
start with a swarm of ~ 0.1-10 km-sized planetesimals. These bodies rapidly grow to ~ 100 - 
1000 km sizes within several Myr in the runaway growth phase (IWetherill and Stewartlll993l ). 
Once icy bodies grow to > 1000 km sizes, they stir the leftover planetesimals to much higher 
velocities. These higher velocities reduce the gravitational focusing factor of the growing 
planets, and the planets enter the much slower oligarchic growth phase. The higher velocities 
also cause more energetic collisions between planetesimals, causing more fragmentation and 
increasing the dust production rate. Becaus e small dust is needed to make debris disks visible 
in the mid-IR, iKenyon and Bromleyl (120081 ) ide ntify the onset of oliga rchic growth with the 
peak in debris disk luminosity. Therefore, the iKenyon and Bromleyl models associate the 
systematically weaker excesses for ^ 5 Myr-old stars with an early stage in planet formation 
where icy bodies are < 1000 km in size and actively growing. Similarly, the models associate 
the stronger excesses surrounding 10-30 Myr-old stars with later stages in planet formation 
where some icy bodies are > 1000 km in size and are growing at slower rates. 



5.3. Constraints on the Primordial-to-Debris Disk Transition for Solar and 

Subsolar-Mass Stars 

The large number of stars with disks in an intermediate stage between primordial disks 
and debris disks allows us to investigate the primordial-to-debris disk transition in detail. 
Our results challenge conventional ideas about two aspects of the primordial-to-debris disk 
transition: the morphology of disks in such an intermediate stage and the timescale for 
transition. 



In agreement with previous Spitzer work, particularly iLada et al.l (120061 ) . our results 
indicate the existence of at least two evolutionary paths from primordial disks to debris disks. 
Many MlPS-detected cluster members have transition disk morphologies with SEDs showing 
weak, optically-thin emission through ^ 5-8 fim but strong, optically-thick emission by ~ 
24 fim. However, many more stars are surrounded by homologously depleted disks whose 
mass is likely lost at all disk radii simultaneously. While 24 fim emission for transition disks 
is stronger than for homologously depleted disks, homologously depleted disks typically have 
stronger < 8 fim emission. It is unlikely that either of these morphologies is the evolutionary 
precursor of the other. Thus, strictly speaking, there is more than one type of 'transition 
disk'. Both morphologies can be thought of as 'evolved primordial disks' since the disk 
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characteristics suggest they are more likely to be composed of first generation dust, not 
second generation dust (i.e. debris disks). 

Our results also place new constraints on the timescale for the transition from primor- 
dial disks to debris disks. According to studies of disks in young, < 3 Myr-old clusters, 
most notably Taurus-Auriga, the paucity of transiti on disks relative to primordial disks 



implies that the transition timescale is 5-. 10 years (ISkrutskie et al.l 119901 : iHartigan et al. 



1990; 



Simon and Prato 



Haisch. Lada. and Lada 



2001 



10; 



995l: Wolk and Waited 19961). Our work and other recent work (e.g. 



Hernandez et al.ll2007al ) indicates that the lifetime of primor- 



dial disks surrounding many young stars is ~ 3-5 Myr. In this picture, the lifetime of the 
primordial disk phase is > 30-50 x the lifetime of the transition disk phase. However, if disks 
spend a factor of 30-50 x longer in the primordial disk phase, then primordial disks should 
always be far more numerous than transition disks. If, on the other hand, the transition 
disk and primordial disk lifetimes are comparable, then transition disks should be equal to 
or greater in number than primordial disks by ~ 5 Myr. 

Based on these data, the traditio nal arguments favoring a 10 ^ yr transition timescale 



are incorrect. For the MlPS-detected Dahm and Hil 



Irwin et al 



enbrandl (120071) sa mple, transition disks 



(120081 ) sample of candidate 



are equal in number to primordial disks. When the 
members is included, transition disks are about twice as numerous as primordial disks. 
Many disks are also homologously depleted and represent a separate evolutionary path from 
primordial disks to debris disks. The co nibined population of M lPS-detected homologously 
depleted disks and transition disks in the iDahm and Hillenbrand! sample are more numerous 
than primordial disks by a factor of ~ 4. When the llrwin et al.l (120081 ) data are included. 



they dominate over the primordial disk population by 6 to 1. Even if the homologously 
depleted disks with the strongest mid-IR emission (IDs 36, 85, 204, 214, and 251) were all 
reclassified as primordial disks, primordial disks would still comprise less than half of the 
total disk population for late-type stars. Because many homologously depleted disks may 
have 24 /im fiuxes lying just below the MIPS detection limit, the population of disks in 
between primordial and debris disks may be even more numerous. 

Our results strongly favor a longer, ~ 1 Myr timescale for the primordial-to-debris disk 
transition for disks with ages ~ 5 Myr. Besides this work, there is evidence in favor of such 
a trend from other recent Spitzer surveys. Spitzer observations of some < 5 Myr-old clusters 
reveal a substantial population of transition disks and homologously depleted disks consistent 
with ~ 1 Myr timescale for the primordial-to-debris disk transition. MlPS-detected disks 
sh owing evidence fo r substantial evolution in 2-3 Myr-old IC 348, la beled anemic disks 



m 



200 



Lada et al.l (120061 ) . comprise ~ 40-45 % of all MlPS-detected disks (ICurrie and Kenyon 
a fraction much higher than that for Taurus-Auriga. SED analysis of these sources 
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shows that about two thirds are homologously depleted disks and one third are transition 
disks (ILada et aLll2006l : ICurrie and Kenyonll2009l ) . Spitzer photometry and IRS spectroscopy 
of low-mass stars in the ~ 1 Myr-old Coronet Cluster also reveal a population of transition 
disks and homolo^ou sly depleted disks whose h igh frequency is in conflict with a < 10^ year 
transition timescale (jSicilia-Aguilar et al.ll2008l ). 



These evolved primordial disks also dominate over primordial disks in at least two 
samples of > 5 Myr-old stars. Most sources with optically-thick mid-IR emission from the 
FEPS sample (ages ~ 3-30 Myr) have SEDs indic ating that their disks' in ner regions are 
cleared of dust to distances probed by ~ 3 .6-8 fim ( iHillenbrand et al.ll2008l ). In four of the 
six stars shown in IHillenbrand et al.l (120081 ) . the SEDs rise from 8-13 nm to ~ 24 /xm as in 



transition disks. The remaining two stars have SEDs showing excess emission weaker than 
that of typical primordial disks and indicative of homologously depleted disks studied here. 
The 2MASS/IRAC survey of h and x Persei (13-14 Myr) fail to detect any stars with clear 
K-band e xcess emission indica tive of primordial disks but find that ~ 4-8% have 8 fim excess 
emission (jCurrie et al.ll2007al ). Wh ile most of these sta rs likely harbor warm debris disks. 



several sources analyzed further in ICurrie et al.l (l2007d ) and one analyzed in ICurrie et al. 



(l2008al ). have accretion signatures and/or mid-IR colors suggestive of evolved primordial 
disks. 



6. Summary and Discussion 

In this paper, we analyze Spitzer data for the ~ 5 Myr-old cluster NGC 2362 to explore 
the primordial-to-debris disk transition and the formation timescale for gas giant planets. 
From the list of confirmed/candidate cluster members, we analyze mid-IR colors, compare 
source SEDs to disk models representating a range of evolutionary states, and compare SEDs 
to signatures of gas accretion to constrain the evolutionary state of disks. Then we compare 
the properties of disks surrounding high-to-low mass stars to constrain the timescale for gas 
giant planet formation. From the sample of early-type, high/intermediate- mass stars with 
debris disks, we investigate the evolution of debris emission from planet formation. Finally, 
we study the morphology of disks evolving from primordial disks to debris disks and place 
limits on the timescale for this transition. 



6.1. Overview of Results 

Our study yields the following major results: 
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• For all stellar masses, NGC 2362 contains very few stars surrounded by primordial disks 
with strong near-to mid-IR emission characteristic of 1 Myr-old stars in Taurus. The 
disk population is dominated by evolved primordial disks, whose SEDs are indicative 
either of so-called transition disks (inner holes/gaps) or homologously depleted disks. 



NGC 2362 lacks any stars with masses > 1.4 Mq that have primordial disks which 
comprise the building blocks of giant planets. These stars also lack evidence for having 
transition or homologously depleted disks. Because most disks surrounding 1 Myr-old 
stars show evidence for protoplanetary disks, gas giant planets surrounding > 1.4 
stars likely form in ^ 1-5 Myr. 

Most disks surrounding solar /subsolar- mass stars are either transition disks or homol- 
ogously depleted disks and thus are likely in the process of dispersing. Some disks may 
be debris disks, indicating that any gas giant planet formation has ceased. Thus, gas 
giant planets may form by ~ 4-5 Myr for stars with a wide range of masses. 

Both canonical transition disks and homologously depleted disks represent an evo - 



lutionary pathway from primordial disks to debris disks (see also iLada et al.l l2006l ) . 
Furthermore, there are a large number of these evolved primordial disks compared to 
primordial disks. This work and previous work suggests that primordial disks have a 
typical lifetime of ~ 3-5 Myr. Thus, the timescale for the primordial-to-debris disk 
transition must be >> 0.1 Myr and may be roughly comparable to the primordial disk 
lifetime (e.g. ~ 1 Myr). A ^ 1 Myr timescale for this evolution is inconsistent with 
published models predicting a rapid transition. 



6.2. On the Primordial-to-Debris Disk Transition 



Despite sample incompleteness in MIPS for the faintest cluster stars, the presence of 
debris disks around the high/intermediate- mass stars and presence of at least many primor- 
dial disks and evolved primordial disks around lower-mass stars in NGC 2362 is consistent 
with the stellar-mass depen dent evolution of primordial disks found in other recent Spitzer 



studies. iLada et al.l (120061 ) showed that disks in 2-3 Myr-old IC 348 exhibit a spectral 



type/stellar- mass dependent frequency of w arm dust from disks with the frequency peaking 
for early M stars. [Hernandez et al.l (j2007al ) also find a stellar-mass dependent frequency of 
warm dust in 3 Myr-old a Orionis. 

Morever, our work suggesting that debris disks may emerge faste r around high/iii t ermed iate- 



mass stars supp o rts th e study of 5 Myr-old Upper Scorpius from ICarpenter et al.l (12006! ). 



Carpenter et al.l (120061 ) find that high/intermediate- mass stars in Upper Scorpius not only 
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have a lower frequency of dis k emission compared to low-mass stars but also have a lower 
16 fim disk luminosity. While ICarpenter et al.l (120061 ) do not model the disk emission from 
high/intermediate- mass stars to demonstrate that these stars are surrounded by debris disks, 
they do note that the weak 16 fim emission and lack of accretion sign atures for these 



stars are consistent with debris disk characteristics. This work and that of ICarpenter et al. 



combined then provide strong evidence that debris disks dominate the disk population of 
high/intermediate- mass stars by ~ 5 Myr. Analysis of the disk po pulation around high-mas s 
stars in IC 348 indicate that debris disks may emerge even sooner (jCurrie and Kenyonll2009l ). 



A 1 Myr timescale for the primordial-to-debris disk transition undermines the vi- 
ability of UV photoevaporation as the primary physical pro cess by which rnost pr imor- 
dial disks disappear as presented in p ublished models from [Alexander et al.l (120061) and 



Clar: 



ce. Gendrin. and Sotomayor (j2001 ). In standard photoevaporation models 



e.g. 



Alexander et al 



20061 ) . once stellar accretion rates drop below the photoevaporative mass loss rate active disk 
dispersal can commence. Once this condition i s met, primordial disks are completely dis- 
persed within ~ 0.01-0.1 Myr. [Alexander et al.l (120061 ) showed that disks dispersing by UV 
photoevaporation follow clear evolutionary tracks in their mid-IR colors: stars with colors 
consistent with bare photospheres and colors consistent with primordial disks should always 
be more frequent than those with intermediate colors. In contrast, our results find that 
disks in an intermediate stage dominate over primordial disks. The inside-out dispersal of 
disks in photoevaporation models is also not consistent with the morphology of homologously 
depleted disks, which comprise the plurality of MlPS-detected disks in NGC 2362. 

The only way in which the NGC 2362 data is consistent with a very rapid transition 
timescale expected from standard photoevaporation models is if the spread in primordial 
disk lifetimes is very narrow (<< 1 Myr) and with a mean value of ~ 4 Myr (R. Alexander, 
pvt. comm.). However, if this were the case for all clusters, plots of the dis k frequency 
vs. time for clusters (e.g. iHaisch. Lada. and Ladall200ll : iHernandez et al.ll2007al ) would look 
more like a step function than a steady decline. Analysis of Spitzer data for other clus- 
ters both younger (Cha I, I C 348, a Ori) and older (TW Hya Association, 25 Ori) suggest 
a larger spread in lifetimes (iLuhman et al.l l2008l : iLada et al.l l2006l : IHernandez et al.l l2007al : 
Low et al.l[2005l : [Hernandez et al.ll2007bh . For each cluster, especially those older than ~ 2—3 
Myr, there exists a wide diversity of disk properties. Some stars around 2-3 M yr old stars 



have very weakly emitting, highly evolved disks or lack any evidence for a disk (ILada et al. 
20061 ) . Conversely, ~ 5-10 Myr old clusters have some st ars with primordial disks, homoL 



ogousl y depleted disks, or canonical transition disks (e.g. ILow et al.ll2005l : IHernandez et al. 



2007bl ). For photoevaporation models to be consistent with observed disk properties, initial 
parameters (e.g. EUV flux, a viscosity parameter) must be fine-tuned to yield a far longer 
disk transition timescale. It is unclear whether these models can modified in such a way. 
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Because our results do not find evidence for a rapid transition timescale, the supposedly 
rapid timescale for all disks as inferred from Taurus data must have alternate explanations. 
Stars in Taurus show an apparently rapid transition timescale because of two main factors: 
1) their youth and 2) binarity. First, previou s Spitzer studies indicate that the typica l 
primordial disk lifetime is at least ~ 3 Myr (ILada et al.l l2006l : [Hernandez et al.l l2007al ). 
which is older than the median age of Taurus (~ 1-2 Myr). Therefore, most disks in Taurus 
(1-2 Myr) are too young to have begun leaving the primordial disk phase. The frequency of 
transition or homologously depleted disks compared to primordial disks for Taurus should 
then be low as observed. While it is possible that some stars (e.g. those in Taurus) make 
the primordial-to-debris transition in ~ 10^ years, our results indicate that a more typical 
transition timescale is far longer. 



Second, binarity affects primordial disk evolution. [Ireland and Kraud (120081 ) show that 
CoKu Tau/4, long regarded a prototypical transition disk, is actually a short-period binary 
system. Its lack of warm dust emission is due to tidal interactions with the faint binary 
companion and not from physical processes typically associated with transition disks (e.g. 
planet formation, photoevaporation). Recent high-contrast imaging of Taurus stars without 
disks indicate that a large number of them also have short-period binaries (A. Kraus et al. in 
prep.). Because binary companions are responsible for clearing these disks, the true number 
of stars lacking a disk because of standard disk dispersal processes (grain growth, accretion, 
planet formation, photoevaporation) is likely far lower than previously thought. Therefore, 
the small fraction of transitional disks in Taurus and other young star-forming regions does 
not necessarily mean that the typical transition timescale due to disk processes is rapid. 



6.3. Implications for the Formation of Gas Giant Planets and Future Work 

Combined with results from previous work, the most straightforward interpretation of 
our results is that gas giant planets have to form sometime between 1 Myr and 4-5 Myr 
around most stars or else they will not form at all. We acknowledge that the strength of 
this conclusion is limited by using indirect tracers for the presence of gas-rich primordial 
disks for many stars. It is also possible that the absolute age calibration of NGC 2362 
is uncertain by several Myr. Nevertheless, a short, 1-5 Myr formation timescale for gas 
giant planets is in strong agreement with timescales inferred from properties of our solar 
system's planets and satellites . For instance, a 1-5 Myr timescale is consistent with the 



results of ICastillo-Rogez et al.l (120071 ) who argue that Saturn formed as early as 2.5 Myr 
after the formation of calcium-aluminum-rich inclusions (CAIs), whose formation may have 
been contemporaneous with the birth of the Sun and the solar nebula. 
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A short formation gas giant formation timescale presents a challenge for core accre- 
tion models of gas gian t planet formation. Many detailed models of planet formation 



[e.g. 



Alibert et al.l |2005| ) assume that circumstellar disks with initial masses comparable 
to the Minimum Mass Solar Nebula are depleted by only factors of ~ 2-5 by ~ 5 Myr. 
Our analysis suggests that p rimordial disks may not last this long. Other models (e.g. 
Dodson- Robinson et al.l l2008l ) claiming to form gas giants within ~ 3-5 Myr rely on sim- 
ple treatments of core growth that neglect important effects such as the evolution of the 
planetesimal veloci ty distribution and co l lision al fragmentation of planetesimals. To date, 
only the model of iKenyon and Bromleyl (120091 ) has successfully modeled the formation of 
gas giant planet cores in < 5 Myr while properly incorporating these first-order effects. The 
relative timescales for c ore growth an d the dispersal of primordial disks may explain trends 
in exoplanet properties (ICurriell2009l ). 



In spite of the failure of most core accretion models to form gas giant planets in less than 
5 Myr while a ccounting for importan t effects, radial velocity surveys indicate that gas giants 



are common (jCumming et al. 



20081). Moreover, g as giants may be most frequent around 



high/intermediate- mass stars (1 Johnson et al.l 120071 ). Our results imply that the formation 



timescales for gas giants around these stars are probably much shorter compared to those 
for lower-mass stars. If gas giants form by core accretion then their formation, especially 
around high/intermediate-mass stars, must be very rapid and efficient. 

Recent and upcoming Spitzer surveys provide an opportunity to explore the primordial- 
to-debris disk transition and planet formation timescales further. Unpublished Spitzer IRS 
data as well as Cycle 5 IRAC and MIPS observations of massive clusters (e.g. Upper Scorpius 
and h and x Persei) strengthen constraints on the primordial-to-debris disk transition at ~ 5 
Myr and probe this transition at later, ~ 10-15 Myr ages. Combined with this study, these 
surveys should result in strong constraints on the evolution of disks from the primordial 
disk-dominated epoch (< 5 Myr) to the debris disk-dominated epoch (> 5-15 Myr) and 
yield the most definitive statements on the timescale for planet formation from the Spitzer 
Space Telescope mission. 



We thank Richard Alexander for useful comments about UV photoevaporation and 
Adam Kraus for very informative discussions on the role of binarity in disk evolution. We also 
thank Carol Grady for pointing out constraints on formation of Saturn from the planetary 
science community. Support for this work was provided by NASA through the Spitzer Space 
Telescope Fellowship Program (TPR). 
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Table 1. Photometry Gatalog for All Stars on the NGG 2362 Field 

Number RA DEC J H Ks [3.6] ct{[3.6]) [4.5] <t([4.51) [5.8] <t([5.8]) [8] ct([8]) [24] ct([24]) 

1 110.0939 -2,-). 1644 13.729 13.486 13.379 13.267 0.011 0.000 0.000 13.18,'5 0.070 0.000 0.000 99.000 0.000 

2 110.1280 -2,"). 3220 12.178 11.34,"i 11.118 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 10.021 0.184 

3 110.0985 -25.1656 12.275 12.206 12.014 11.964 0.005 0.000 0.000 0.000 0.000 0.000 0.000 99.000 0.000 

4 110.0921 -25.1472 16.153 15.419 15.431 16.186 0.039 0.000 0.000 0.000 0.000 0.000 0.000 99.000 0.000 
6 110.0960 -25.1629 14.618 14.142 14.064 13.977 0.017 0.000 0.000 13.434 0.086 0.000 0.000 99.000 0.000 



Note. — The table includes photometry for all stars on the NGC 2362 field with a detection in IRAC and/or MIPS. Values of 99.000 or 0.00 in 
the photometry column and 0.00 in the photometric uncertainty columns denote sources that were not observed in a given filter. 
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Table 2. Photometry Catalog for NGC 2362 Cluster Members 

ID RA DEC V R I J H K,, [3.6] <t([3.6]) [4.5] ct([4.51) [5.8] o-([5.8]) [8] o-([8]) [24] ct([24]) 

1 109.5817 -24.9550 19.430 18.470 16.800 15.220 14.620 14.370 14.130 0.050 14.110 0.070 14.125 0.148 14.054 0.138 99.000 0.000 

2 109.5905 -24.9578 17.420 16.450 15.650 14.670 13.950 13.770 13.580 0.040 13.590 0.050 13.491 0.086 13.705 0.097 99.000 0.000 

3 109.5909 -24.8615 18.170 17.170 15.950 14.600 13.860 13.580 13.240 0.030 13.060 0.040 12.745 0.047 12.089 0.022 7.649 0.009 

4 109.5942 -25.0283 19.680 18.530 16.920 15.540 14.910 14.520 14.280 0.050 14.250 0.070 13.979 0.129 13.885 0.117 99.000 0.000 

5 109.6002 -24.8997 18.890 17.840 16.420 14.930 14.410 14.010 14.000 0.060 13.850 0.060 13.819 0.113 13.992 0.130 99.000 0.000 

Note. — Tho table includes photometry for NGC 2362 cluster members identified bv lDahm and HillenbrandH2007ll . Values of 99.000 or 0.00 i 
the photometry column and 0.00 in the photometric uncertainty columns denote sources that were not observed in a given filter. 



Table 3. Photometry Catalog for NGC 2362 Candidate Members from llrwin et al.l (120081 ) 



ID RA DEC VI ,7 H Ks [3.6] ct([3.61) [4.5] ct([4.51) [5.8] ct([5.8]) [8] <t([8]) [24] o-([24]) 

1 109.5003 -25.1379 20.515 17.525 15.985 15.346 14.938 14.835 0.030 14.831 0.041 14.773 0.265 14.637 0.273 99.000 0.000 

2 109.5145 -25.1769 16.786 15.526 14.564 13.971 13.854 0.000 0.000 13.793 0.020 0.000 0.000 14.005 0.135 99.000 0.000 

3 109.5160 -25.1295 21.548 18.795 99.000 99.000 99.000 16.067 0.079 16.044 0.111 0.000 0.000 0.000 0.000 99.000 0.000 

4 109.5212 -25.1615 20.936 18.074 16.399 15.728 99.000 15.392 0.045 15.228 0.056 14.925 0.308 14.879 0.369 99.000 0.000 

5 109.5234 -25.2406 22.386 18.925 16.617 15.570 14.916 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 9.800 0.064 

Note. — The tabic includes photometry for candi date NGC 2362 cluster members identified bv llrwin et al.l i2008>) . Values of 99.000 or 0.00 : 
the photometry column and 0.00 in the photometric uncertainty columns denote sources that were not observed in a given filter. 



Table 4. Dust Removal Timescales and Dust Masses for Early- Type Stars with MIPS 24 

fim Excesses 

ID J-Ks Spectral Type L*/Lq M/Mq Ttj^st (AU) tpr (yr) s,„ai Disk max(Mdust) 

(Estimated) (Estimated) (Estimated) (200, 120 K) (10, 30 AU) (/^m) Model (xIO^^OMq) 



31 0.11 A5 


30.0 


2.3 


10.7,29.8 


2.3x10-1, 2.1x10^ 


14.8 


3020154-5-9 


14 


144 0.45 G5 


10.0 


2.0 


6.2,17.2 


7.0x10'!, 6.3x10^ 


5.8 


3008807-5-9 


1.1 


303 0.01 AO 


35.3 


2.5 


11.6,32.3 


2.0x10'', 1.8x10^ 


16.0 


3017274-5-9 


2.2 


314 -0.05 B9 


56.5 


2.7 


14.7,40.9 


1.2x10'*, 1.1x10^ 


23.8 


3017274-5-9 


2.2 



Note. — Constraints on the dust removal timescales and dust masses for early-type sources. The Poynting-Robcrtso n drag 
timcscale is calculated assuming 10 /jm-sized grains. Disk models listed represent the file number from the iRobitaille et al.l 
l l2006i1 grid; dust masses for those models are listed as M^^gt. 
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Table 5. Evolutionary State of Confirmed/ Candidate NGC 2362 Members with Evidence 

for Disk Emission 



ID ST A„ Range EW(H„) log(Mjji3fc/M0) Ri„_„i„ (AU) /3 Disk Type Notes 



Dahm and Hillenbrand (2007) sample with MIPS detections 



3 


M2 


0-0.65 


NA 


-4, 


,51± 0.48 


4, 


,10 


1 


.05 


Transition 




36 


K5 


0-0.53 


-6.7 


-3, 


,02 


± 1, 


.06 


0, 


,07 


1 


.02 


Homologously 


Depleted 


41 


Kl 


0-0.33 


-0.1 


-7, 


,39 


± 0, 


.42 


0, 


,12 


1 


.05 


Homologously 


Depleted 


63 


MO^ 


0.31? 


-2.9 


-6, 


,62 


± 1, 


.02 


0, 


,06 


1 


.06 


Homologously 


Depleted 


85 


M2 


0-0.93 


-45.4 


-3, 


,57 


± 0, 


.93 


0, 


,02 


1 


.14 


Homologously 


Depleted 


111 


K7 


0-0.45 


-18.8 


-2, 


,07 


± 0, 


.36 


0, 


,08 


1 


.13 


Primordial 




139 


M2 


0-0.8 


NA 


-2, 


,09 


± 0, 


.39 


0, 


,11 


1 


.13 


Primordial 




168 


K3 


0.55-0.9 


-0.2 


-5, 


,63 


± 


.85 


5 


,20 


1 


.09 


Transition 




177 


M2 


0-0.2 


-3.0 


-3, 


12 


± 


.97 


0, 


,06 


1 


.11 


Homologously 


Depleted 


187 


M2 


0.03-1.3 


-117.7 


-3, 


,58 


± 


.7 


0, 


,27 


1 


.04 


Primordial^ 




194 


K5 


1.02-1.8 


-27.5 


-3, 


,91 


± 


.57 


7, 


,57 


1 


.12 


Transition 




196 


MO^ 


0.31? 


NA 


-5. 


,95 


± 


.54 


3, 


,63 


1 


.04 


Transition 




202 


M3 


0-1.08 


-2.6 


-2 


41 


± 


.12 


0, 


,03 


1 


.18 


Primordial 




204 


M3 


0-1.05 


-10.4 


-3, 


,87 


± 


.85 


0, 


,03 


1 


.16 


Homologously 


Depleted 


213 


M3 


0-0.8 


-10.2 


-4, 


,35 


± 1 


.03 


0, 


,04 


1 


.10 


Homologously 


Depleted 


214 


K7 


0-0.43 


-3.5 


-3, 


,63 


± 1 


.00 


0, 


.06 


1 


.10 


Homologously 


Depleted 


219 


MO 


0.78-1.68 


-5.9 


-4, 


24 


± 1 


.14 


0, 


.03 


1 


.07 


Homologously 


Depleted 


229 


Ml. 5 


0-0.95 


-3.8 


-3, 


,81 


± 1 


.05 


0, 


,03 


1 


.07 


Homologously 


Depleted 


251 


M2^ 


0.31? 


NA 


-3, 


,98 


± 


.82 


0, 


,02 


1 


.15 


Homologously 


Depleted 


267 


M2^ 


0.31? 


NA 


-4, 


,35 


± 


.79 


0, 


,02 


1 


.12 


Homlogously Depleted 


295 


M2^ 


0.31? 


NA 


-4, 


,36 


± 1 


.12 


1 


,29 


1 


.08 


Transition 










Irwin et 


al. (2008) sample 


with MIPS 


detections 




IR-1S5 


M2^ 


0.31? 


NA 


-4, 


12 


± 1, 


.05 


0, 


,75 


1 


.07 


Transition 




IR-459 


M2^ 


0.31? 


NA 


-4, 


,59 


± 1, 


.57 


0, 


,05 


1 


.06 


Homologously 


Depleted 


IR-526 


M2^ 


0.31? 


NA 


-3, 


,62 


± 0, 


.76 


0, 


,03 


1 


.14 


Primordial 




IR-583 


M2l 


0.31? 


NA 


-3, 


77 


± 0, 


.94 


0, 


,6 


1 


.12 


Transition 




IR-663 


M2l 


0.31? 


NA 


-4, 


,41 


± 0, 


.93 


0, 


,02 


1 


.10 


Homologously 


Depleted 


IR-732 


M2^ 


0.31? 


NA 


-3, 


,30 


± 1, 


.27 


0, 


,04 


1 


.05 


Homologously 


Depleted 


IR-750 


M2^ 


0.31? 


NA 


-3, 


10 


± 0, 


.72 


4, 


,32 


1 


.09 


Transition 




IR-809 


M2^ 


0.31? 


NA 


-4, 


,93 


± 1, 


.07 


0, 


,03 


1 


.08 


Homologously 


Depleted 


IR-817 


M2^ 


0.31? 


NA 


-3, 


22 


± 0, 


.63 


1 


,95 


1 


.11 


Transition 




IR-851 


M2^ 


0.31? 


NA 


-3, 


,25 


± 0, 


.51 


0, 


,21 


1 


.14 


Transition 




IR-931 


M2^ 


0.31? 


NA 


-4, 


24 


± 1, 


.27 


0, 


,39 


1 


.03 


Transition 




IR-958 


M2^ 


0.31? 


NA 


-5, 


,09 


± 1, 


.19 


0, 


,04 


1 


.07 


Homologously 


Depleted 


IR-1046 


M2^ 


0.31? 


NA 


-2, 


,9 ± 0.95 


0, 


,04 


1 


.13 


Primordial 




IR-1091 


M2^ 


0.31? 


NA 


-4, 


,39 ± 1.64 


1 


,45 


1 


.06 


Transition 






Disk Candidates from 


Dahm 


and Hillenbrand (2007) 


without MIPS detections 


12 


M3 


0.9-1.72 


-20.2 


















Homologously Depleted 3 


27 


M3 


0.2-1.03 


-7.2 


















Homologously Depleted 3 


39 


M2^ 


0.31? 


NA 


















Homologously Depleted 3 


100 


K5 


0.3-1.1 


NA 


















Primordial 


4 


227 


M2^ 


0.31? 


< 


















Homologously Depleted 3 



Note. — Footnote (1) identifies stars whose spectral types are estimated from optical/near-infrared colors. Footnote (2): IDs 187 and IR-526 
have low disk masses but are classified as primordial based on their mid-IR fluxes (compared to the median Taurus SED). Note (3) in the last 
column — Stars with probable 2—3 cr detections at MIPS 24 (j.m. Note (4) in the last column — MIPS detection is blended with another star. The 
range in extinction, A^, is calculated by comparing the observed V-I colors to intrinsic colors from [Kenyon and Hartman 

3 119951 1. assuming a 1 

subclass uncertainty in the spectral type. Negative values for EW(Hq,) identify stars with Hq, emission. Sources lacking a measured EW{Ha,) have 
NA listed in the EW{Ha) column. The column log(M^^sfc /Mq ) refers to the unweighted average and standard deviation of all models that provide 
a good fit; R^^i min is the minimum inner radius of all models with good fits. 




Fig. 1. — Histogram of IRAC detections with 2MASS counterparts. The 5.8 fxm (dash-three 
dots) and 8 //m (sohd) source counts peak at m[IRAC] ~ 15 and 14.25, respectively. The 
3.6 iim (dash-dot) and 4.5 iim (dashed hne) both peak at m[IRAC] ~ 15.5. 
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Fig. 2. — (top) Histogram of MIPS 24 fj,m detections (solid line) and sources with 2MASS 
counterparts (filled region), (bottom) Error distribution of 24 //m detections with 2MASS 
counterparts. Sources with photometric errors slightly greater than the main distribution 
typically he closer to the image edges. 
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Fig. 3. — False-color image of NGC 2362 made with the PBCD mosaic images at 4.5 jjm 
(blue), 8 fim (green), and 24 fim (red). The bright star r CMa is shown near the image 
center. 



Fis- 4. — J-K,s vs. Ks-[5.8 ] (left) and K5-[8] color-color diagrams for sources in the 



Dahm and Hillenbrand! (120071 ) membership list using photometry derived in this paper. 
Sources lacking MIPS detections are shown as triangles. Sources with MIPS detections 
are shown as dots. Overplotted are the locii of photospheric colors (dotted lines) from 
the Kurucz-Lejeune stellar atmosphere models as calculated by the SENS-PET tool. Most 
sources with IRAC excesses have MIPS detections. 




Fig. 5. — Sa me as Figure HI except fo r the li st of candidate members from llrwin et al.l (120081 ). 
As with the iDahm and Hillenbrandl (120071 ) sample, most sources with IRAC excesses have 
MIPS detections. 
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Fig. 6. — J vs. Ks-[24] (left) and J-K ,, vs. K..-[24] (right) co l or-col or diagrams for con- 
firme d /candidate cluste r members from iDahm and Hillenbrand! (120071 . top) and candidates 
from Irwin et al. ( 20081 . bottom). The leftward arrows identify Kg- [24] upper limits for 
sources lacking MIPS detections. 
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-[24] 



Fig. 7. — The K<j-[5.8,8] vs. Ks-[24] color-color diagrams for the iDahm and Hillenbrand 



( 2007 ) sample (top). Overplotted are the colors from the median Taurus SED (jPurlan et al. 
2006h reddened by E(B-V) =0.1 (square), (top-left) The K,-[5.8] vs. K,-[24] color-color 



diagram. Sources with 24 fim upper limits are shown as grey left-pointing arrows. MIPS- 
detected sources (black dots) are shown overp lotted with la photome tric errors, (top-right) 
The Kg- [8] vs. Ks-[24] color-color diagram for iDahm and Hillenbrand! sources. The symbols 
are the same as in the le ft-hand panel. (bottom)Ks- [5.8,8] vs. Ks-[24] color-color diagrams 
for the llrwin et al.l (120081 ) sample. 
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Fig. 8. — SEDs of early /intermediate- type sta rs with masses ^1.4 Mq. Overp lotted are the 
best-fitting SEDs from iRobitaille et al.l (l2006l . solid line) along with the stellar photosphere 
(dotted line). Also listed is the J2000 source position (in degrees). Photometric uncertainties 
are smaller than the symbol sizes. Inferred dust masses are 10"^° Mq to 10"^ Mq, consistent 
with debris disks. 
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Fig. 9. — Representative spectral energy distributions of late-type NGC 2362 members and 
candidate members with MIPS 24 fim excesses illustrating the range of disk evolutionary 
states. Overplotted are the median Taurus SED (dotted line) with upper and lower quartiles 
(dot-dashed line). A terrestrial zone debris disk model is shown as a dashed line with mid-IR 
flux slightly greater than the photosphere (solid line, connected by open squares). The MIPS 
5(7 detection limit is shown as a horizontal grey dashed line. The derived disk evolutionary 
states are (clockwise from the top left) primordial, homologously depleted, transition, and 
debris. 



1 10 
A (,um) 



1 10 100 1000 

Disk inner radius (R l^) 




-4.35± 0.79 



2000 5000 10' 2x10' 

Stellar temperature (K) 



10-" 




10 ' 10 ^ 10 

Disk mass (Mg) 



Fig. 10. — SED analysis of ID-267 using the lRobitaille et al.l (120061 ) models to fit the observed 
optical/infrared flux (top- left), constrain the inner disk radius and flaring index /? (top-right), 
constrain the stellar temperature (bottom left), and constrain the disk mass (bottom rig ht). 
The solid line in the top-left plots shows the best-fitting SEDs from the iRobitaille et al.l grid 
with other possible models overplotted as faint, grey lines. In the plot of /? vs. Rj„ (top-right 
panel), black dots denote the flaring index and inner disk size (in units of the sublimation 
radius) for the models shown in the top-left panel. The grey dots show values for the entire 
grid. The black curves in the bottom panels define the mean and standard deviation of 
the stellar temperature and disk masses from the best-fitting models, the latter assuming 
a gas-to-dust ratio of 100. The hatched regions show the actual disk masses for each SED 
model shown the top-left panel. The distributions of temperatures and disk masses for the 
entire grid are shown in grey. 
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Fig. 11. — Spectral energy distribution s for primordial/evolved primordial disk candidates 
from the iDahm and Hillenbrandl (120071 ) sample lacking MIPS detections. Symbols are the 
same as before. 
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Fig. 12. — Amplitude of 24 /xm excess around BAF stars vs. time for 4-40 Myr old clusters. 
Besides NGC 2362 (filled dots), we include Orion Obla and b (10 and 5 Myr; triangles), rj 
Cha (8 Myr; stars), the /? Pic Moving Group (12 Myr; squares), h and x Persei (13 Myr; 
diamonds), Sco-Cen (5, 16, and 17 Myr; asterisks), NGC 2 232 (25 My r ; inve rted triangles), 
and NGC 2547 (38 Myr; small crosses). HR 4796A, from ILow etaP koO^ is also shown 
as a cross at t~ 8 Myr and [24]-[24]p/iot ~ 5. The data are consistent with a rise in debris 
emission from 5 My r to 10 Myr, a peak from ~ 10 Myr to 20 Myr as shown previously by 
Currie et al.l (l2008al ). For stars older than 40 Myr, th e distribution of 24 fim luminosity 
defines an envelope which is consistent with a t~^ decline ( iRieke et al.ll2005l ). The dashed line 
shows the mean excess. Predictions for the decline in emission from a steady-state collisional 
evolution model (dash-three dots) follow a t~^ decline. 
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A. A Comparison between our IRAC photometry and that from Dahm and 

Hillenbrand (2007) 

Our primary motivation for redoing IRAC photometry o f NGC 2362 i s to id entify candi- 
date cluster stars with mid-IR excesses at < 10 fim from the llrwin et al.l (120081 ) catalog and 
to construct full SEDs from V band through 24 /im to constrain the origin of the emission 
in all cluster s tars. To test the accura c y of o ur IRAC p hotometry we compare i t with pho- 
tometry from iDahm and Hillenbrand! (120071 ) . Because iDahm and Hillenbrand (120071 ) only 
publish photometry for cluste r members, our analysis in this section is limited to the 337 
Dahm and Hillenbrand! (l2007l ) (candidate) members. 



Figures [13] and [H] contrast our photometry with iDahm and Hillenbrand! (12007! ). In 
Figure [131 we show the difference in magnitude for each IRAC channel as a function of 
(our) magnit ude. Photometry from both g roups appears to have little/no zero-point offset. 
However, the lDahm and Hillenbrand! (!2007l ) photometry exhibits a deflection towards fainter 
magnitudes for m(IRAC) > 13 in all filters relative to our photometry. This effect is especially 
pronounced in the [5.8] and [8], which clearly show a non-gaussian scatter. 

Figure UM shows that the mid-IR colors also have systematic differences. The top panels 
compare the distributions of J-K^ vs. Ks-[8] colors. We have far fewer sources with with 
unphysical K5-[8] colors (< -0.2): the blue limit of the K.,- [8] distribution tracks the photo- 
spheric locus (solid line) very well. . The Ks-[8] colors from iDahm and Hillenbrand! (120071 ) 
deflect blueward of the photospheric locus by up to 0. 5-1. Our photometry include s only 
two stars with Ks-[8] ~ -0.2 compared to seventeen in iDahm and Hillenbrand! (120071 ). The 
[5.8] vs. [3.6]-[5.8] distribution for our sources (middle panels) also contains fewer stars with 
very blue ([3.6]-[5.8] < -0.5) colors and exhibits a smaller dispersion in color, especially for 
stars fainter than m [5.8] = 13. 



In the bottom panels of Figure [HJ we quantify the differences in photometry at 4.5 jjm 
and 8 /im by fitting a line (dash- dot) to the 8 fim to 4.5 fim flux ratio (log( F[8]/F[4.5])) as a 



function of J-K^ color (see also [Carpenter et al. 



2006!: 



Carpenter et al. J2006h and 



Currie et al. 



Currie et al. 



2008br) for so urces with 



J2008bh . the best- 



5cr detections in both filters. As in 
fit line is determined by a linear least squares fit with 4cr clipping to remove outlying sources, 
incl uding stars with infrared exces s. The dispersion and rms residuals are then recalculated. 
The iDahm and Hillenbrandl (12007! ) photometry has a 3a dispersion in log(F[8]/F[4.5]) of ~ 
0.23 and residuals in the linear fit of ~ 5.5%. Five stars with larger photometric uncertainties 
have negative fiux ratios that lie outside the dispersion (dashed lines). Our photometry has 
a disp ersion of 0.16 and 4.3% fit res iduals with one outlying star. Furthermore, four stars 
in the iDahm and Hillenbrandl (12007! ) photometry have extremely negative fiux ratios, lying 
5-6 a fainter than the photospheric locus. If the distribution in fiux ratios followed a normal 
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distribution, we would expect, at most, one such source. Our photometry contains just one 
3(7 outher with a negative flux ratio. 

The differences in photom etry are primarily due to the s maller aperture radius employed 
here (2 pixels vs. 3 pixels in iDahm and Hillenbrand! 120071 ). The sky background may be 
overestimated from the larger sky annulus needed for a 3 pixel aperture in densely populated 
regions, which results in artificially faint sources whose colors deflect blueward. Differences 
in image processing methods may also explain differences in photometry. Differences in 
the aperture radius and image processing methods have led to significant discrepancies in 
photometry for other Spitzer cluster data. In these case s, more accurate ph otometry was 
also obtained using a smaller aperture radius (e.g. Cha I; iLuhman et al.ll2008l ). 
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Fig. 1 3. — Difference between our IRAC photom etry and that from iDahm and Hillenbrand 
( 2007 ) as a function of our IRAC magnitude. The Dahm and Hillenbrandl ( '20071 ) photometry 
exhibits a deflection towards fainter magnitudes as a function of survey depth (5 (IRAC) < 
0). 



Fig. 14. — IRAC colors from iDahm and Hillenbrandl (l2007l . left columns) and this work 
(right c olumns), (top) The J-K. vs. Kg- [8] color-color diagram for sources in common 
between iDahm and Hillenbrandl (120071 ) and this work, (middle) [5.8] vs. [3.6]-[5.8] color- 
magnitude diagrams, (bottom) The 8 /xm to 4.5 fim flux-ratio diagram for the same samples. 
The best-fit line (dash-dot line) and 3 a dispersion above and below the best-fit line (dashed 
lines) are overplotted. The lDahm and Hillenbrandl (120071 ) photometry yields many stars with 
unphysically blue 2MASS/IRAC colors (K,-[8] < -0.2; [3.6]- [5.8] < -0.5) and unphysically 
small flux ratios (< -1.15). The photometry from this work lacks many blue colors and 
exhibits less far scatter, especially at fainter magnitudes in the 8 fim channel. 
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B. Atlas of SEDs of MlPS-detected (candidate) cluster members 
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Fig. 15. — Atlas of spectral energy distributions of late-type NGC 2362 members with MIPS 
24 fim excesses. Overplotted are the median Taurus SED (dotted line) with upper and 
lower quartiles (dot-dashed line). A terrestrial zone debris disk model is shown as a dashed 
line with mid-IR flux slightly greater than the photosphere (solid line, connected by open 
squares) . The MIPS 5a detection hmit is shown as a horizontal grey dashed hne. 



-57- 











S'tronb IFJaC' ' ' 10^.698 


1 0"'* 








m2 —24.9760 — 










1 77 : 












10"'^ 
















































\x 


10"™ 










10"=' 





















10"'* 




Weak'iRA'C ' ' ' 


109.706 - 






-^i^-.^ k5 


-24,9404 : 








194 : 






• 




10"'^ 












• 
















X 




10"™ 








10"=' 



















S'trong IFJaC' ' ' 


109.710 : 

-24.9655 - 






202 




xv • 


• 




















1 10 
X (microns) 




S'trong lf?AC' ' ' 


109.72' : 
-24.9437 '- 






215 




\x • 








• - 










10 




^ m2 


109.707 ■ 
-24,9740 
195 ~- 




\X 






\X 










1 10 
X (r^icrons) 




S'trong II^AC' ' ' 


109.712 '- 
-24,8930 : 
204 




^C-K • "V--.-- 

Xx 






\. X 











X (microns) 





s'trong If^AC' ' ' 




109.722 




"Hfe ^ k7 




-24.9656 - 








214 : 




X. •""X.;----, y""" 
































\x 












• - 




x^^x 























Fig. 16. — Same as Figure 
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Fig. 1 8. — Spectral energy distributions for MlPS-detected candidates from the llrwin et al, 
(120081 ) ■ Symbols are the same as in the preceeding figures. 



